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PREFACE 


This  report  is  taken  from  the  thesis  of  M.  W.  Beckstead 
(Ph.D. ,  University  of  Utah,  1965).  The  work  was  performed 
under  sponsorship  of  the  Air  Force  Office  of  Scientific  Re¬ 
search,  Grant  AF-AFOSR  446-63.  It  was  presented  at  the  AIAA 
Third  Aerospace  Sciences  Meeting  in  New  York  City,  January, 
1965. 


ABSTRACT 


Non-acoustic  combustion  instability  has  been  examined  in  an 
uncatalyzed,  a  catalyzed,  and  two  aluminized  composite  propellants. 
These  propellants  were  studied,  burning  cigarette  fashion,  in  a  burner 
capable  of  operating  at  values  of  L*  as  small  as  5  cm*  It  vaa  °t" 
served  that  the  frequency  of  the  pressure  oscillations  varied  with 
the  value  of  L*,  frequency  decreasing  with  increasing  L* .  The  data 
were  correlated  by  plotting  frequency  versus  the  reciprocal  of  V, 
yield*. .g  a  series  of  constant  pressure  curves.  This  pressure  effect 
was  eliminated  by  using  dimensionless  variables,  allowing  all  of  the 
data  for  a  given  propellant  to  be  correlated  along  the  same  line. 

A  one-dimensional  model  is  proposed  that  considers  sinusoidal 
perturbations,  allowing  for  growth  of  the  disturbance.  The  pressure, 
the  burning  rate,  the  distributed  temperature  in  the  propellant,  and 
the  energy  flux  from  the  burning  gases  are  the  quantities  perturbed. 
The  analysis  is  developed  utilizing  an  energy  balance  on  the  solid 
phase  that  accounts  for  energy  accumulation  in  the  solid  and  a  mass 
balance  on  the  rocket  chamber  that  considers  mass  accumulation  in  the 
gas  phase.  The  assumption  is  made  that  the  energy  flux  from  the  gas 
phase  either  oscillates  in  phase  with  the  pressure  or  is  slightly 
Out  of  phase. 

The  analysis,  thus  developed,  reduces  to  an  expression  relating 
the  reciprocal  of  the  dimensionless  L*  to  the  dimensionless  frequency 

iv 


in  terms  of  the  growth  factor  and  propellant  parameters .  The  principal 
parameters  involved  are  the  activation  energy,  the  mean  surface 
temperature,  and  the  heat  of  gasification  for  the  propellant.  With 
the  proper  choice  of  these  parameters  the  theoretical  analysis  yields 
results  that  agree  quantitatively  with  the  observed  experimental  results. 
The  theoretical  results  also  appear  to  be  in  qualitative  agreement 
with  moat  of  the  experimental  and  theoretical  results  obtained  by 
other  investigators. 

The  restate  of  this  investigation  do  not  answer  all  of  the  many, 
unanswered  questions  in  non-acoustic  instability  but  can  be  considered 
significant  in  three  respects.  First,  an  experimental,  L*-frequency 
dependency  was  observed  in  non-acoustic  instability,  and  was  corre¬ 
lated  through  the  use  of  a  dimensionless  frequency  and  L*.  Secondly, 
a  mathematical  expression  has  been  derived  that  agrees  quantitatively 
with  the  observed,  experimental  results,  subject  to  certain  assumptions. 
Finally,  the  results  of  this  investigation  can  be  used  to  correlate 
qualitatively  some  of  the  apparently  unrelated  observations  of  other 
investigators . 
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CHAPTER  I 


INTRODUCTION 

Combustion  Instability  In  solid  propellants  has  received  con¬ 
siderable  attention  over  the  yeaurs  because  of  the  destructive  effects 
It  has  on  a  rocket  motor.  It  may  cause  only  a  Blight  vibration  or 
change  In  the  thrust  of  the  rocket ;  It  may  result  In  complete  des¬ 
truction  of  the  rocket  motor;  or  In  some  cases  It  may  extinguish  the 
burning  propellant  and  thus  terminate.  None  of  the  manifestations 
of  combustion  Instability  In  a  rocket  motor  is  desirable. 

Combustion  instability  as  applied  to  solid  propellants  can 
roughly  be  divided  into  two  categories:  acoustic  and  non-acoustic. 

TOie  wave  equation,  written  for  the  combustion  chamber  of  a  rocket 
motor,  will  give  the  acoustic  frequency  of  the  chamber,  and  pressure 
oscillations  corresponding  to  this  frequency  can  be  considered  as 
acoustic  oscillations  or  acoustic  instability.  The  majority  of  the 
published  literature  concerning  combustion  instability  in  solid  pro¬ 
pellants  falls  Into  this  category.  A  theoretical  analysis  of  acoustic 
instability  advanced  by  Hart  and  McClure  [17]  has  been  quite  extensively 
developed  and  appears  at  the  present  to  be  the  most  readily  accepted 
theory  in  that  area. 

Oscillations  have  also  been  reported  [see  for  example,  1,  3; 21, 26, 

33,  35,  39  ]  which  have  occurred  at  frequencies  much  lower  than  the  acoustic 
frequency  of  the  chamber  in  which  they  have  appeared.  This  type  of 
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Instability  has  been  termed  non-acoustic  or  low  frequency  instability 
(typically  0-500  cps),  and  is  that  which  will  be  considered  in  the 
present  investigation  (the  designation  "non-acoustic  instability" 
being  abbreviated  to  NAl).  This  phenomenon  will  not  be  referred  to 
as  low  frequency  instability  because  the  author  feels  that  this  termi¬ 
nology  can  be  deceiving  and  is  not  as  descriptive  as  that  of  non-acouatic 
instability.  The  sureas  of  interest  in  NAI  can  be  classified  as  chuffing 
and  non-acoustic  oscillatory  combustion.  Chuffing  is  the  periodic 
ignition  and  extinguishment  of  propellant  (see  Figure  l),  whereas  oscil¬ 
latory  instability  is  characterized  by  periodic  or  quasi-periodic 
pressure  disturbances  about  a  mean  pressure  (see  Figure  1  also). 

It  has  been  stated  earlier  that  there  is  not  a  great  deal  of 
published  literature  in  the  field  of  non-acoustic  instability,  and 
most  of  what  is  available  has  appeared  in  very  recent  years.  There 
are  several  reasons  why  more  research  is  being  done  on  NAI  now  than 
before.  One  reason  is  that  metal  additives  have  been  incorporated  into 
propellants  in  an  effort  to  suppress  acoustic  instability.  A  degree  of 
success  has  been  attained,  but  the  presence  of  metal  additives  has 
apparently  aggravated  the  mechanism  that  initiates  NAI.  A  second 
reason  appears  to  be  due  to  the  fact  that  many  of  the  newer  and  more 
modem  rockets  operate  at  relatively  low  pressures  compared  with  the 
pressures  employed  in  esurlier  rockets,  and  these  lower  pressures  seem 
to  encourage  the  appearance  of  NAI.  Still  another  reason  stems  from 
the  fact  that  motors  have  increased  in  size  to  the  extent  that  the 
characteristic  acoustic  frequencies  have  correspondingly  been  reduced 
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Figure  1.  Non-Acoustic  Combustion  Instability. 
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until  they  now  are  of  the  same  order  of  magnitude  as  non-acoustic 
oscillations  in  smaller  motors.  The  problem  thus  arises  of  possible 
cjupling  between  acoustic  and  nor -acoustic  instability. 

Some  of  the  earliest  reported  work  wherein  one  investigator  recog¬ 
nized  the  possibility  of  a  different  mechanism  than  that  which  exiBts 
in  acoustic  instability,  was  that  of  Huffington  in  1951*  [21  ]  and  1956 
[10].  He  observed  and  reported  on  the  phenomenon  of  chuffing  using 
cordite  as  propellant.  Experimental  results  indicated  that  the  frequency 
of  oscillations  increased  with  pressure.  He  also  observed  that  the  dis¬ 
tance  the  surface  regressed  during  a  chuff  was  constant  for  a  given 
pressure.  The  Frank-Kamenetskii  thermal  explosion  theory  tl3 1 
used  to  explain  the  experimental  results.  This  theory  is  based  on  the 
assumption  that  exothermic  reactions  are  occurring  throughout  a  homo¬ 
geneous  medium,  causing  self -heating.  For  a  critical  thickness  of 
material,  the  reactions  will  reach  a  run-away  condition  and  an  explosion 
will  occur.  In  Huffington' s  work  the  distance  the  propellant  regressed 
during  a  chuff  was  proposed  as  being  the  critical  thickness  and  data 
were  correlated  on  this  basis.  While  some  double-base  propellants 
apparently  do  have  exothermic  reactions  occurring  throughout  the  solid 
phase,  many  propellants  for  which  the  existence  of  such  reactions  is 
questionable  have  exhibited  chuffing  and  non-acoustic  instability  [26]. 
This  fact  has  discouraged  extensive  use  of  the  thermal  explosion  theory. 

However,  Angelus  and  Young  [3,35]  using  modified  double-base 
propellants  have  followed  a  course  indicated  by  the  thermal  explosion 
theory,  by  correlating  the  induction  time  between  chuffs  with  the 
average  pressure  during  a  chuff.  Angelus  was  one  of  the  first 
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frequencies  within  which  the  non-acoustic  oscillations  lie  at  a  given 
pressure  and  for  a  specific  propellant.  Data  taken  with  a  large 
experimental  burner  having  acoustic  modes  which  lie  in  the  non-acoustic 
frequency  range  indicate  that  low  frequency  acoustic  oscillations  fall 
approximately  within  the  preferred  band  of  frequencies  defined  by  non¬ 
acoustic  oscillations.  This  means  that  the  coupling  between  acoustic 
and  non-acoustic  instability  which  was  mentioned  earlier  can  actually 
take  place  and  could  be  a  very  serious  problem  with  large  solid  pro¬ 
pellant  boosters . 

A  theoretical  model  for  low  frequency  instability  in  solid  pro¬ 
pellant  rocket  motors  has  been  proposed  by  Akiba  and  Tanno  [1].  They 
assumed  steady-state  conditions  and  then  through  use  of  standard 
perturbation  techniques  examined  the  limits  of  combustion  stability. 

The  equations  which  they  used  were  the  unsteady  heat  conduction 
equation 


ST 

3t  = 


d2T 

a  — 5  +  r 
3x2 


St 

S 


(l.i) 


(for  a  definition  of  all  symbols,  see  the  Table  of  Nomenclature  in 
Appendix  A),  a  surface  energy  balance 


k  H  +  -  k  H  -  prQ  (1.2) 

g  s+  6  8 

(where  Hs  is  the  temperature  gradient  at  the  surface,  and  Q  is  the 
heat  evolved  by  surface  reaction)  and  an  equation  relating  the  depen¬ 
dence  of  the  temperature  gradient  at  the  surface  (on  the  gas  phase 
side)  to  the  pressure 
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Is  v  -  (Tf  -  *.)  k6  f  (1-3) 

S 

/  ML*  \ 

where  K  is  a  chamber  time  constant  (K  «  —  )• 

g  g  CDRTf 

By  substituting  a  perturbed  temperature,  pressure,  burning  rate,  and 
temperature  gradient  at  the  surface,  assuming  that  the  gaseous  phase 
responds  much  faster  than  the  perturbations  under  consideration  (low 
frequency  oscillations),  and  eliminating  the  steady-state  terms  in  the 
customary  manner,  they  obtained  a  chamber  transfer  function  and  a  com¬ 
bustion  transfer  function.  The  Nyquist  stability  criterion  was  then 
applied  to  these  transfer  functions  in  order  to  determine  stability 
limits  for  the  burning  propellant.  The  theoretical  results  Indicated 
that  the  reciprocal  of  the  motor  free  volume  (l/L*)  was  proportional 
to  pressure  raised  to  twice  the  burning  rate  exponent  (n).  The  experi¬ 
mental  results  which  were  presented  for  a  double-base  propellant 
appeared  to  agree  quite  well  with  the  theoretical  predictions.  However, 
it  1b  difficult  to  determine  the  precise  meaning  of  their  results 
because  no  explanation  was  made  of  experimental  apparatus  or  the 
criteria  used  to  determine  the  recorded  data  points .  One  drawback  to 
the  theoretical  results  of  this  method  of  determining  a  stability  limit 
is  that  it  predicts  a  critical  frequency  at  the  stability  boundary 
only  and  does  not  give  any  relationship  for  the  dependency  of  the 
frequency  of  the  oscillations  for  other  conditions,  nor  was  any 
reference  in  the  paper  made  as  to  the  agreement  between  this  predicated, 
critical  frequency  and  the  actual,  experimental  frequencies . 
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Sehgal  and  Strand  [35]  have  modified  the  model  proposed  by  Akiba 
and  Tanno  by  using  a  mass  balance  on  the  rocket  in  place  of  Equation 
(1.3)  in  order  to  obtain  a  chamber  transfer  function.  The  mass  balance 
which  was  used  was 

St  <V=h>-  -  CDPchAt  (l-k) 

where  is  a  nozzle  discharge  coefficient.  The  rest  of  the  analysis 
followed  the  same  outline  as  that  of  Akiba  and  Tanno  with  the  exception 
that  the  paper  in  general  was  more  complete.  The  experimental  data 
were  obtained  from  regressive  burning  grains  which  burned  until  the 
low  pressure  extinction  limit  was  reached  for  the  existing  conditions . 
The  remaining  amount  of  propellant  was  measured  and  the  value  of  L* 
for  the  chamber  at  this  point  was  then  calculated.  The  data  deter¬ 
mined  in  this  manner  are  not  strictly  compatible  with  the  theoretical 
assumptions  because  the  analysis  concerned  itself  with  instability 
triggered  by  small  perturbations,  whereas  the  data  presented  [see 
also  2}  corresponded  to  flame  extinction  after  a  period  of  finite 
pressure  oscillations .  Anderson,  Strehlow  and  Strand  [2 ]  have 
reported  additional  data  of  this  type  for  various  polyurethane-based 
propellant  compositions.  They  observed  that  as  the  percentage  of 
aluminum  in  the  propellant  increased,  so  did  the  "slope  of  the  data 
although  the  burning  rate  exponent  remained  approximately  the  same. 

Thus  for  increasing  aluminum  content  there  appeals  to  be  an  increasing 
discrepancy  in  the  analysis .  Studies  were  also  carried  out  which 
indicated  that  incomplete  combustion  was  not  a  serious  problem  for  the 
range  in  which  data  were  taken. 
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Prom  the  above  discussion  of  the  research  which  has  been  carried 
out  in  the  area  of  NAI,  it  can  be  seen  that  different  reports  are 
sometimes  contradictory  and  that  none  are  conclusive .  Although  a 
considerable  effort  has  been  made  in  very  recent  years  and  is  continuing 
at  the  present,  neither  a  precise  mechanism  nor  a  completely  satis¬ 
factory  mathematical  model  has  been  prepared  to  explain  the  occurrence 
of  non-acoustic  instability.  The  present  investigation  does  not  pur¬ 
port  to  answer  the  questions  but  the  intent  is  that  it  will  add  to  the 
available  knowledge  and  understanding  of  the  subject  in  a  positive  way. 

At  the  time  when  the  present  study  was  begun  (1962),  very  little 
information  concerning  NAI  was  available  in  the  literature.  Most  of 
the  above  work  had  Just  been  stented  or  was  in  progress  but  had  not  yet 
been  reported.  Therefore,  the  initial  goal  that  was  established  was  to 
learn  more  about  the  general  nature  of  the  problem.  As  experimental 
procedures  were  developed  and  results  began  to  be  positive  in  nature, 
this  goal  was  more  precisely  defined  as  determining  the  stability 
limits  for  a  particular  group  of  propellants.  An  uncatalyzed,  a 
catalyzed,  and  an  aluminized  composite  propellant  were  chosen  in  an 
attempt  to  obtain  fairly  representative  data  without  studying  a  wide 
variety  of  propellants.  These  propellants  have  been  designated 
as  Utah  0,  F,  and  TF  propellants  respectively,  and  their  exact 
compositions  and  physical  properties  are  given  in  Appendix  B. 

Experiments  were  carried  out  with  the  hope  that  from  the  observations, 
a  mathematical  approach  considering  the  propellant  properties  and  the 
chamber  conditions  might  be  established  which  would  allow  the  prediction 
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of  a  stability  boundary  for  the  onset  of  instability.  In  an  attempt 
to  fulfill  this  goal  a  mathematical  model  is  proposed  and  developed 
in  Chapter  II. 

Negative  initial  experimental  results  from  shaped  grains  led  to 
the  decision  to  employ  a  burner  capable  of  operating  at  very  small 
free  volume  (low  L*).  Subsequently  an  end  burner  similar  to  one  used 
at  the  Naval  Ordnance  Test  Section  [12]  was  designed  and  constructed. 
The  data  presented  in  the  following  chapters  were  all  obtained  using 
this  burner  in  Its  various  configurations. 


CHAPTER  II 


THEORETICAL.  ANALYSIS 

This  chapter  proposes  a  model  of  non-acoustic  Instability, 
states  the  attendant  assumptions,  and  presents  the  derivations  of  the 
resultant  equations,  but  before  proceeding  it  seems  appropriate  to 
comment  briefly  on  the  general  nature  of  the  system  involved. 

A  composite  solid  propellant  consists  of  a  crystalline  oxidizer 
sind  a  polymeric  binder  which  acts  as  a  fuel.  The  oxidizer  is  often 
added  as  a  bimodal  mixture  in  order  to  improve  the  mixing  qualities 
of  the  propellant  while  it  is  being  compounded.  The  resultant  material 
is  a  heterogeneous  aggregate  of  Intimately  mixed  crystals  of  odd  shapes 
and  sizes  and  the  fuel  binder.  This  heterogeneous  nature  of  solid  pro¬ 
pellants  makes  It  difficult  if  not  impossible  to  describe  the  problem 
of  combustion  instability  mathematically  without  making  many,  rather 
drastic  assumptions. 

A.  A  Model  of  Non-Acoustic  Combustion  Instability 

Several  assumptions  relating  to  the  propellant  and  the  rocket 
chamber  have  been  made  by  others  considering  similar  problems  concerned 
with  propellant  combustion.  These  assumptions  will,  therefore,  be 
considered  as  "usual  simplifying  assumptions"  and  will  not  be  discussed 
in  detail  here.  These  include  assuming  that:  the  solid  is  mathematically 
semi-infinite,  the  physical  properties  of  the  propellant  and  the  burned 
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gas  do  not  change  with  temperature  over  the  range  of  interest,  the 
solid  is  homogeneous  and  isotropic,  the  gas  obeys  the  ideal  gas  law, 
isentropic  flow  exists  through  a  sonic  nozzle,  and  burning  occurs 
uniformly  on  the  assumed  flat  surface  of  the  propellant.  The  assump¬ 
tion  of  homogeneity  is  probably  the  most  questionable  of  the  assumptions 
and  will  be  considered  more  in  detail  at  a  later  time.  Other  assump¬ 
tions  peculiar  to  the  problem  at  hand  will  be  pointed  out  as  they  are 
made. 

Other  investigators  [18,25]  have  made  order  of  magnitude  cal¬ 
culations  which  indicate  that  the  mechanism  of  energy  transport  in 
the  solid  phase  must  be  considered  in  an  analysis  of  low  frequency 
combustion  instability  (which  is  directly  applicable  to  the  problem 
being  considered).  The  following  paragraphs  include  a  brief  discussion 
of  the  physical  considerations  involved  in  analyzing  the  energy  distri¬ 
bution  in  the  propellant  and  its  response  to  low  frequency  perturbations. 

An  examination  of  the  steady-state  temperature  distribution  in  a 
burning  solid  propellant  appears  to  be  a  good  starting  point  for  the 
following  discussion.  A  schematic  sketch  of  this  postulated,  idealized 
temperature  profile  has  been  included  as  Figure  2.  Temperature  is 
taken  as  the  ordinate,  and  the  the  distance  into  the  propellant  is 
measured  in  the  positive  x-direction.  The  origin  is  considered  as 
being  stationary  and  at  the  burning  surface  with  the  propellant  moving 
toward  it  at  a  velocity  equal  to  the  linear  burning  rate.  The  area 
beneath  the  profile  is  proportional  to  the  sensible  energy  of  the  pro¬ 
pellant.  It  is  this  energy  and  its  response  to  external  (gas-phase) 


X*0  X  increasing 


Figure  2.  Idealized  steady-state  temperature  profile  in  a 
burning  solid  propellant. 
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perturbations  which  will  be  considered  in  the  following  analysis . 
Consider  the  steady-state  heat  conduction  equation 


a 


0 


Applying  the  boundary  conditions 


at  x  -  0,  T  *  i 

O 


(2.1) 


at  x  *  »,  T  =  Tq 

gives  the  solution  to  (2.l)  as  [for  example,  see  38,  p.  79,  or  35) 

r 

_  —  x 

T  -  T0  +  (Ts  -  TQ)  e  °  (2.2) 

Therefore,  the  steady-state  temperature  profile  is  dependent  upon  the 
average  surface  temperature,  the  thermal  diffusivity,  the  burning  rate, 
umi  the  initial  uniform  temperature,  TQ.  It  is  generally  accepted 
that  the  surface  temperature  varies  with  pressure  (increasing  as  pres¬ 
sure  increases),  though  the  explanation  is  disputed  [see  22,4,23). 
However,  the  burning  rate  dependence  upon  pressure  is  much  greater 
thon  that  of  the  surface  temperature,  though  related  to  it,  so  that 
the  burning  rate  is  the  more  important  psLrameter  in  determining  the 
shape  of  the  profile,  and,  therefore,  the  amount  of  energy  stored  in 
the  solid. 

At  steady-state  conditions  the  burning  rate  of  the  propellant 
is  regulated  by  the  rate  of  energy  feedback  from  the  burning  gas.  In 
the  following  analysis  this  energy  will  be  considered  as  an  external 
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flux  and  will  be  equated  to  the  rate  of  erergy  consumption  by  the 
propellant  due  to  heating  and  reactions  at  the  surface. 

The  manner  In  which  this  energy  transport  is  accomplished  will 
be  considered  as  being  mainly  due  to  thermal  conduction  with  radiation 
effects  being  of  a  secondary  importance  [6,37]-  It  has  been  pointed 
out  [6,7]  that  the  effects  of  radiation  should  become  greater  at  very 
low  pressures  and  low  burning  rates  and  therefore*  these  conditions 
will  be  considered  further  at  a  later  time  and  are  discussed  in  some 
detail  in  Appendix  G. 

If  the  system  thus  described  is  subjected  to  a  small  perturbation, 
then  the  equilibrium  between  these  energy  fluxes  will  be  disturbed  as 
will  the  burning  rate,  the  mass  of  gas  generated,  the  mass  flux  through 
the  nozzle,  and  the  chamber  pressure.  In  responding  to  a  perturbation 
of  this  nature  the  solid  is  able  to  act  as  a  thermal  capacitance  in 
the  sense  that  a  surplus  of  energy  can  be  stored  momentarily  and 
released  later,  causing  the  rate  of  energy  consumption  to  be  out  of 
phase  with  the  rate  of  energy  supply.  If  the  energy  fluxes  and  the  above- 
mentioned  mass  fluxes  react  to  each  other  in  such  a  manner  that  the 
perturbation  grows,  then  the  system  is  unstable  for  the  given  conditions. 

The  thermal  capacitive  effect  mentioned  above  can  be  readily 
associated  with  the  time-dependent  difference  between  the  energy  supply 
and  consumption  in  the  solid.  This  idea  is  developed  rigorously  in 
the  following  analysis.  The  relationship  between  this  thermal  capaci¬ 
tance  and  the  other  variables  of  the  system  that  could  describe  instability 
is  not  readily  obvious .  However,  a  time-dependent  mass  balance  describing 
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the  mass  flow  through  the  rocket  chamber  will  be  utilized  to  provide 
a  relationship  of  this  nature. 

It  is  hoped  that  this  brief  discussion  of  the  phenomenological 
occurrences  taking  place  in  a  rocket  motor  during  unstable  combustion 
will  enable  the  reader  to  better  understand  the  mathematical  develop¬ 
ments  of  the  following  sections.  In  the  analysis  that  follows,  a 
standard  perturbation  technique  will  be  applied  wherein  the  derived 
equations  will  be  subjected  to  small  perturbations.  The  first-order 
perturbation  terms  will  be  kept  after  steady-state  terms  are  subtracted 
out,  and  second-order  terms  will  be  neglected  altogether. 

B.  A  Perturbed  Energy  Balance 

Other  investigators  [1,15,35,36]  considering  the  problem  of 
combustion  instability  in  solid  propellants  have  utilized  a  surface 
energy  balance  in  their  analyses.  These  developments  have  considered 
the  temperature  gradient  at  the  surface  as  being  indicative  of  the 
energy  flux  to  the  surface,  and  have  applied  Fourier's  law  of  heat 
conduction  to  the  burning  surface  of  the  propellant  (for  example, 
see  Equation  (l.2)).  In  the  present  analysis  a  slightly  more  general 
approach  will  be  taken  wherein  all  of  the  thermal  energy  stored  in  the 
solid  will  be  considered,  and  the  feed-back  flux  from  the  gas  phase 
will  not  be  tied  strictly  to  the  surface  temperature  gradient. 

The  over-all  energy  balance  for  the  propellant  from  the  sur¬ 
face  to  a  station  deep  in  the  propellant  is 

<f  +  %  '  =,  *0>  -  »p  r  <’  *  'p  T„>  *  It  I0  DP  %  TdX  (2'3) 
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or,  by  simplifying 

f  ■>  Pp  r  fq  +  Cp  (Ts  -  T0)]  +  CpPp  f”  |f  dx  (2.3a) 

where  f  is  the  energy  flux  to  the  surface  of  the  propellant,  and  q  is 
the  net  energy  absorbed,  per  unit  mass,  by  the  gasification  process 
and  could  possibly  be  either  positive  or  negative. 

It  is  assumed  that  the  deviations  from  the  steady  state  are  of  a 
sinusoidal  nature.  Figure  3  shows  the  growth  of  pressure  oscillations 
on  a  chuff,  and  it  can  be  seen  that  they  are  periodic  and  approximately 
sinusoidal.  It  will  therefore  be  assumed  that  the  perturbations 
of  all  of  the  quantities  are  sinusoidal.  The  temperature,  burning 
rate,  flux,  and  pressure  can  now  be  written  as  the  sum  of  their  steady- 
state  value  plus  a  small  perturbation  that  has  a  sinusoidal  time 


dependence.  The  equations  are 

T(x)  =  T(x)  +  Tg  6(x)  e^i+T|^CDt  (2.4a) 

r  -  r  [1  +  3  e(l+Tl)a3t]  (2. kb) 

f  =  f  [1  +  ♦  e^1+T|^Dt]  (2.4c) 

P  =  F  [1  +  «  e(i+Tl)cDt]  (2.4d) 

The  exponent  from  the  time-dependent  portion  of  (2.4)  can  also  be 
expressed  in  the  following  manner. 

(i  +  t])  cot  =  (a  +  iu>)  t  (2.5) 


where  the  term  on  the  right-hand  side  of  the  equation  can  be  considered 
as  a  complex  frequency  or  preferably  as  the  frequency  plus  the  growth 
factor  (ag). 


-  18  - 


Run  No.  42-063 

L*  i  23  cm. 

Ordinate  scale  -  10  psi/rrajor  division 
Abscissa  scale  =  0.025  sec/major  diviscn 


Figure  3.  Periodic  nature  of  non-acoustic  oscillations 
superimposed  on  a  chuff . 
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It  should  be  noted  here  that  the  pressure,  the  flux,  and  the 
temperature  sure  not  necessarily  in  phase  with  the  burning  rate. 
Consequently,  three  of  the  four  quantities  p,  9,  ♦,  and  Jt  are  in  general 
complex.  In  this  analysis,  p  is  taken  as  read  and  the  others  are  complex. 

Substitution  of  Equations  (2.4)  into  Equation  (2.3a),  subtracting 
out  the  steady-state  terms,  and  discarding  the  second-order  perturbation 
terms  yields 


♦ 

P 


1  + 


c  T 


+ 


C  CD  T 
2 _ 8 

rH 


(n  +  l) 


dx 


(2.6) 


This  equation  establishes  a  relationship  between  the  ratio  of  the  flux 
perturbation  to  the  burning  rate  perturbation  a is  modified  by  energy 
accumulation  in  the  solid.  In  order  to  evaluate  the  integral  in  (2.6) 
an  equation  relating  the  temperature  to  the  burning  rate  as  a  function 
of  distance  is  needed.  The  unsteady  heat  conduction  equation  will  be 
employed  to  derive  such  a  relationship. 


C.  The  Perturbed,  Unsteady  Heat  Conduction  Equation 

Energy  transport  in  the  solid  phase  will  now  be  considered. 

Heat  conduction  will  be  regarded  ais  one  dimensional,  bulk  phase  reactions 
will  be  neglected,  and  the  frame  of  reference  will  be  the  same  as  in 
Figure  2.  The  unsteady  heat  conduction  equation  is 


-  St 
+  r  5^ 


(2.7) 


Substituting  Equations  (2.4a)  and  (2.4b)  into  (2.7)  and  regrouping 
results  in 


0  = 


a*2  L  "  ax2 


+  rT 


d0 

dx 


+  ?p  §  -  fs<D0(i+ri) 


(2.8) 
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The  first  group  of  terms  in  (2.8)  is  nothing  more  than  the  steady-state 
heat  conduction  equation  and  is  equal  to  zero  (see  Equation  2.1 ). 
Simplifying  leaves 


d2e  r  de  cp(i-m)  Q  n  £r_  dT  (2.9) 

dx2  adx  ai  ^ 

o 

The  steady-state  temperature  gradient,  which  is  needed  to  solve 
Equation  (2.9),  can  be  obtained  readily  by  differentiating  Equation 
(2.2).  Substituting  this  result  into  (2.9)  produces  the  following 
second-order,  nonhomogeneous  equation 


^0  ,  ?x  46  .  e 

„  2  dx  a 

dx 


(2.10) 


The  boundary  conditions  for  the  system  are 

at  x  =  0,  9(0)  =  0. 

at  x  =  oo,  9(oo)  =  0 


where  0Q  is  the  amplitude  of  surface  temperature  perturbation.  The 
solution  to  Equation  (2.10)  can  be  obtained  by  applying  the  standard 
techniques  for  solving  differential  equations  [for  example,  see 
reference  37,  P*  79) •  The  solution  is 


e(x)  'fo  +  ^s~T0^  1  e-X(l+a+lt)x  _  U(*s  T0^  e  -2Xx 
P  "  7(i+r,)fs  -*  7(i+T!)T8 


(2.11) 


where  X,  a,  a,  and  t  are  defined  as  follows 


X  =  L 

2a 


(2.11a) 
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(2.11b) 
(2.11c) 

(2. lid) 

Both  o  and  r  are  Introduced  to  aid  in  simplifying  the  form  of  the  equation. 
The  dimensionless  frequency,  7,  is  similar  to  the  frequencies  and  time 
constants  used  by  Green  [15  ],  Hart  and  McClure  [l8],  and  others  [1,35]* 

The  motive  for  perturbing  the  unsteady  heat  conduction  equation 
vaa  to  provide  a  relationship  between  the  perturbed  temperature  and 
burning  rate  as  a  function  of  distance.  Equation  (2.1l)  provides  that 
relationship,  and  can  now  be  combined  with  Equation  (2.6)  to  evaluate 
♦,/p.  The  actual  integration  of  the  integral  in  Equation  (2.6)  is 
carried  out  in  Appendix  C  along  with  further  algebraic  manipulations. 

The  results  are  (see  Equation  (C-7)) 


|  -  1  +  0  +  Gx  [  (t](ct+1)  +  t]  +  Y(o+l  -  |)+i  [G1(o+l-qT)  -  YT ] 

(2.12) 

c  T  0-. 
where  G  =  “ 

G  -  7CPTs  fo  _  2L  0 

G1  2HS  3  2S  G 

- 

Y  =  HS 

2  2 

and  S  =  t  +  (o+l) 
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Equation  (2.12)  describes  the  coupling  between  the  energy  flux 
and  the  burning  rate  perturbations,  with  the  capactive  effect  of 
energy  storage  in  the  propellant  built  in.  It  is  now  necessary  to 
obtain  a  relationship  involving  the  motor  parameters;  a  mass  balance 
will  be  utilized  to  accomplish  this. 


D.  A  Perturbed  Mass  Balance 

Many  investigators  [for  example,  see  1,2,26, 28, 33> 35 3  have 
observed  that  the  mean  pressure  and  the  L*  of  a  rocket  chamber  have 
an  effect  on  NAI.  This  section  will  be  devoted  to  deriving  a  relation¬ 
ship  involving  these  variables. 

Sehgal  and  Strand  [35]  used  a  mass  balance  for  the  rocket  chamber 
to  relate  these  variables,  and  the  same  mass  balance  (see  Equation 
(1.4))  will  be  used  here  but  will  be  developed  in  a  slightly  different 
manner.  The  mass  balance  for  a  rocket  motor  containing  burning  pro¬ 
pellant  and  exhausting  through  a  sonic  nozzle  is  [see  5]. 


ft  (oevch' 


»PrAb  - 


V. 

c* 


(2.13) 


where  C*  is  the  "characteristic  velocity"  for  the  products  of  combustion 
as  described  in  texts  on  internal  ballistics  [5 J -  Differentiating  the 
accumulation  term  of  (2.13),  assuming  that  ideal  gas  conditions  exist 
in  the  chamber,  and  recognizing  that  the  time  rate  of  change  of  the 
volume  is  equal  to  the  area  of  the  burning  surface  multiplied  by  the 
burning  rate  gives 


Vch  dp 
RTf  dt 


r  ^  (pp  -  Pg) 


(2.14) 
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Recognizing  that  the  density  of  the  gas  Is  negligible  compared  to  that 
of  the  solid,  dividing  (2.14)  by  At,  and  regrouping  yields 


L*  dP  „  „  A  P 

RT^  dt  “  S  *  Pp  C* 


(2.15) 


Substituting  the  perturbed  pressure  and  burning  rate  (Equations  2.4b 
and  2.4d)  into  (2.15),  and  eliminating  the  steady-state  terms  in  the 
usual  manner  leaves  the  following 


oj  (i+n)  P*  -  pp  Kjj  r  0  - 
Recognizing  from  steady-state  considerations  that 

and  simplifying,  results  in 


where 


£ 


■=  1  +  tj7S  +  i  75 


--  -2 
C*  r  .. 

UoRT~  L 


(2.16) 


(2.17) 


The  expression  0/x  (the  ratio  of  the  burning  rate  perturbation  to 
pressure  perturbation)  is  called  the  response  function.  The  importance 
of  the  response  function  is  quickly  realized  when  one  considers  the 
overall  goal  of  an  investigation  such  as  the  present.  If  the  delicate 
balance  between  the  pressure  and  the  burning  rate  is  disturbed,  the 
question  as  to  whether  conditions  will  return  to  an  equilibrium  state, 
or  run  away,  is  paramount.  The  response  function  appears  to  be 
intimately  involved  in  the  answer  to  this  question. 
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E.  The  Response  Function 

From  the  definitions  of  the  perturbed  quantities  (Equations  (2 A)) 
it  can  be  seen  that  for  a  positive,  negative,  or  zero  growth  factor,  t), 
the  following  conditions  are  realized*  For  t)  >  0,  there  exist  growing 
perturbations  and  unstable  conditions;  for  rj  ■  0,  there  is  neutral 
equilibrium  representing  a  stability  boundary;  and  for  r\  <  0  there 
are  decreetsing  perturbations  and  stable  conditions.  From  Equation 
(2.17)  it  can  be  seen  that  for  the  above  conditions  the  real  part  of 
the  response  function  will  be  greater  than  unity,  equal  to  unity,  and 
less  them  unity,  respectively.  It  is  interesting  to  note  that  Hart 
and  McClure  in  their  analysis  of  the  coustic  instability  problem 
used  the  response  function  as  a  stability  (or  amplification)  criterion 
[17],  and  at  low  frequencies  the  stability  limit  is  defined  as  the 
point  where  the  real  part  of  tie  response  function  1b  equal  to  unity 
[16]. 

If  an  independent  means  of  evaluating  the  growth  factor  were 
available,  there  would  be  no  need  of  further  considerations  in  order 
to  define  a  stability  boundary.  However,  at  the  present  the  only 
direct  method  of  approximating  the  growth  factor  is  from  experimental 
data,  and,  therefore,  further  deliberations  are  necessary. 

F.  Derivation  of  the  Frequency-L*  Dependency 

In  Section  C  an  expression  for  the  ratio  of  the  perturbed  flux 
to  burning  rate  was  derived  (see  Equation  (2.12))  and  in  Section  D  an 
expression  for  the  ratio  of  perturbed  burning  rate  to  pressure  was 
derived  (Equation  (2.17)).  The  present  task  is  to  combine  these 


ti 
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equations  In  a  meaningful  manner  which  can  be  used  to  describe  the 
physical  situation.  As  both  (2.12)  and  (2.17)  contain  the  burning 
rate  perturbation,  the  obvious  is  to  combine  them,  eliminating  the 
burning  rate  effect.  This  can  be  done  by  multiplying  (2.12)  by  (2.17), 
the  result  being 

£  -  (l  +  T]7t)  Re  (^/3)  -  7t  im  (e/p)  +  i  (rfc  Re  (e/p)  +  (l  +  nrt)  Im(e/P)] 

(2.18) 

e 

where  Re  (♦/?)  *  1  +  G  +  [*)(o+l)  +t ]  +  y  (<x+l  -  g) 

and  Im  (#/p)  -  G1  (o+1-t|t)  -  Yt 

From  (2.l8)  It  can  be  seen  that  the  problem  Is  now  reduced  to  that  of 
evaluating  the  ratio  of  the  flux  to  pressure  perturbations.  The  following 
paragraph  considers  a  relationship  of  this  nature. 

It  was  mentioned  in  Chapter  I  that  a  large  experimental  burner 
having  acoustic  modes  that  lie  in  the  non-acoustic  range  has  been 
constructed  and  utilized  quite  extensively  by  the  Naval  Ordnance  Test 
Station.  This  burner  exhibits  a  definite  "preferred  frequency"  regime 
for  given  propellantB  outside  of  which  oscillations  are  not  observed 
[12,19,20,28].  It  appears  reasonable  to  postulate  that  the  solid 
phase  is  supplying  the  mechanism  necessary  to  cause  oscillations  at  a 
non-acoustic  frequency  and  that  the  acoustic  modes  of  the  burner  are 
supplying  the  necessary  coupling  to  allow  these  oscillations  to  be 
maintained  (the  acoustic  mechanism  replacing  the  low  L*  or  non-acoustic 
coupling  mechanism).  Some  work  has  been  done  with  this  burner  in 


-  26  - 


studying  the  phase  relationship  between  the  oscillations  of  the  light 
intensity  irom  the  burning  gases  and  the  pressure  oscillations.  It 
has  been  observed  that  the  pressure  is  not  always  in  phase  with  the 
emitted  light.  This  phase  shift  has  been  observed  in  both  double-base 
propellants  [12]  and  composite  propellants  [29].  During  a  typical  run 
in  this  burner,  oscillations  begin  to  grow  as  the  pressure  increases 
with  the  pressure  oscillations  lagging  those  of  the  light  intensity. 
This  phase  relationship  passes  through  zero  as  the  amplitude  of  the 
oscillations  passes  through  a  maximum.  The  oscillations  die  out  with 
the  pressure  leading  the  light.  From  this  it  appears  that  the  response 
of  the  solid  is  a  maximum  when  the  pressure  and  the  light  intensity 
are  in  phase.  If  the  observed  light  intensity  is  associated  with  the 
energy  flux,  it  would  appear  most  profitable  to  Investigate  small, 
positive  or  negative  phase  shifts,  or  a  zero  phase  shift. 

Equation  (2.l8)  can,  therefore,  be  modified  to  include  the  phase 
angi for  the  flux-to-pressure  perturbation.  Doing  this  and  solving 
for  the  dimensionless  L*  gives  » 


,  _ _ F  tan  £  -  1 _ 

75  =  F  (1  -  Tjtan  £)  +  (q  +  tan  5) 


where  F 


_  Re  (0 


Im  (VP 


(2.19) 


and  tan  £  s  Ref^xj  ^  being  the  phase  angle) 

In  Equations  (2.l8)  and  (2.19),  o,  t,  T  ,  G,  Y,  and  S  are  functions 
of  7,  i),  the  propellant  parameters,  an£  tan  £.  In  the  following 

section  two  methods  of  approximating  0Q/p  from  the  propellant  properties 
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will  be  presented,  and  1;  will  be  assumed  that  tan  ?  will  be  small  or 
can  be  approximated.  It  can  therefore  be  seen  from  Equation  (2.19) 
that  the  frequency  is  a  function  of  L*  and  the  growth  term. 

Although  the  problem  appears  to  have  been  reduced  to  a  solution 
for  n  involving  only  the  two  variables  7  and  £,  there  still  exists 
a  problem  of  a  practical  nature.  Physical  properties  and  burning 
rate  data  for  propellants  can  be  obtained  with  a  certain  amount  of 
accuracy,  but  numerical  values  for  chemical  properties  such  as  the 
activation  energy,  heat  of  gasification  for  the  propellant,  heats  of 
reactions  for  the  surface  reactions,  and  the  surface  temperature  are, 
at  best,  crude  approximations.  Realizing  this,  ranges  of  values  for 
the  above  three  uncertain  parameters  were  obtained  from  the  literature, 
where  possible  (average  values  were  taken  from  References  22  and  2U  )  • 

A  parametric  study  wa is  th°n  carried  out  to  determine  the  best  agree¬ 
ment  of  the  experimental  data  with  the  approximated  parameters  for 
various  values  of  tan  £.  A  discussion  of  these  results  cam  be  found 
in  Chapter  V. 


G.  Evaluation  of  ®Q/P 

Kinetic  Limited  Case.  -  A  relationship  between  the  surface 
temperature  and  the  burning  rate  must  be  obtained  in  order  to  evaluate 
9q/@.  The  following  Arrhenius  expression  is  used 


-e/rt 

r  **  A  e  8 


(2.20) 


Substituting  the  perturbed  burning  rate  and  surface  temperature 
expressions  (Equations  (2.k))  into  (2.20),  assuming  that  the  perturbations 
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are  small  compared  to  unity,  and  solving  for  0Q/p,  one  finds  [see  reference 

351 

RT 

9o/6  .  (2.21) 

Thus,  if  the  value  of  the  surface  temperature  is  determined  by  kinetic 
limitations.  Equation  (2.2l)  provides  an  approximation  for  0Q/p  which 
can  then  be  used  with  Equation  (2.19)  giving  the  dimensionless  L* 
as  a  function  of  the  dimensionless  frequency  and  the  growth  term. 

The  assumption  that  the  surface  temperature  determines  the 
burning  rate  has  been  popular.  It  has  been  utilized  in  most  of  the 
theories  of  unstable  combustion  [1,15,17,35)  with  the  numerical  value 
of  the  activation  energy  ranging  between  16  Kcal/mole  and  50  Kcal/mole 
in  the  various  publications. 

Thermodynamic -Limited  Case. --Recently  it  has  been  proposed  that 
the  surface  temperature  of  a  composite  propellant  burning  at  low 
pressures  iB  controlled  by  the  equilibrium  decomposition  of  solid 
ammonium  perchlorate  [4,23,24].  The  chemical  equation  describing  the 
decomposition  is 

NH^CIO^  -  NH3  +  HC10U  (2.22) 

The  heat  of  sublimation  for  this  reaction,  calculated  from  vapor  pressure 
measurements,  has  been  reported  as  56  Kcal/mole  [22,23) •  Experiments 
carried  out  by  Fowling  and  Smith  [24)  indicate  that  the  value  of  the 
surface  temperature  is  fixed  by  the  pressure  up  to  60  psia  and  probably 
higher.  This  represents  the  lower  range  of  working  pressures  involved 
in  the  present  study,  and  therefore  warrants  further  investigation. 
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Assuming  that  the  surface  temperature  is  determined  by  the  equili¬ 
brium  vaporization  of  Equation  (2.22),  the  Clausius -Clapeyron  equation 
can  then  be  utilized  in  the  following  form 


P  AH 


s 


T 


(2.23) 


s 

Substituting  the  perturbed  burning  rate  and  surface  temperature 
expressions  (Equation  (2.4))  into  2.23),  assuming  that  the  perturbations 
are  small  compared  to  unity,  and  solving  for  0q  one  finds  (see  Appendix 
D,  Equation  (D-2)) 


9 


0 


(2.24) 


B 

The  fact  that  (2.24)  is  a  real  equation  and  not  complex  indicates 
that  eQ  is  in  phase  with  x  which  follows  logically  from  the  assumption 
that  the  surface  temperature  is  fixed  by  the  pressure.  Instead  of 
attempting  to  solve  for  0q/(3,  the  expression  for  0q  (Equation  (2.24)) 
will  be  substituted  directly  into  Equation  (2.12),  resulting  in 


|  -  1  +  Y2  (|)  +  Y1  (|)  h(cr+l)+T]  +  Y  (l+a-|)  +  i  [Y]L  (|)(1+0-tit)  -  Yt] 

(2.25) 

t,  m  2 

y  c  R  T 

where  Y1  5 

s 

c  R  T  2 
Y2  S  PHAH 

s 


Multiplying  this  equation  by  the  response  function  in  order  to  eliminate 
the  burning  rate  perturbation  as  was  done  to  obtain  Equation  (2.l8) 


gives 
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|  =  (1+TirS)  [l+Y  (1  +  o-  §)]  +  Y2  +  Yx  [tj  (o+l)  +  '1  +  YT75 

+  1  J^Sr  [1  +  Y  (l  +  O  -  |  -  tjt)]  +  (l  +a-T]T)  -  Yt  j  (?.26) 
Following  the  same  argument  that  led  to  Equation  (2.19)  yields 

tan  C  {l  +  Y  (l+o-§)  +Yg  [1  +  g§  (t)(o+1)+t)]}  +  Yt  -  Yg  fg  (1+p-t|t) 
7£  l+Y  (l+a  -  tan  £  [yt  +  q [l+Y(l+a-|) j|  ^ 

If  the  assumption  is  made  that  the  surface  temperature  is  in  phase 
with  the  pressure  rather  than  the  burning  rate,  then  Equation  (2.27) 
applies  in  place  of  Equation  (2.19). 

As  to  the  superiority  of  the  kinetic  assumption  over  the  thermo¬ 
dynamic  assumption,  or  vice  versa,  Nachbar  and  Williams  [22]  report 
that  the  evidence  of  their  studies  indicates  that  neither  mechanism 
can  be  verified  as  the  only  mechanism,  nor  eliminated  as  an  incorrect 
mechanism.  Therefore,  both  assumptions  have  been  considered  in  the 
present  analysis  and  results  from  using  each  assumption  in  connection 
with  the  parametric  study  will  be  presented  in  Chapter  V. 

H.  Summary  of  Theoretical  Considerations 

The  results  of  the  analysis  made  in  this  chapter  are  essentially 
contained  in  Equations  (2.19)  and  (2.27).  Both  of  these  equations 
predict  that  the  frequency  of  the  oscillations  is  related  to  the  L*  of 
the  motor,  and  furthermore,  in  dimensionless  form  this  relationship 
is  independent  of  the  pressure.  This  result  is  quite  revolutionary 
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considering  that  although  it  is  quite  generally  accepted  that  low  values 
of  L*  are  characteristic  for  non-acoustic  instability  and  various 
low  frequencies  have  been  observed  experimentally,  to  the  author's 
knowledge  a  correlation  between  L*  and  frequency  has  not  been  reported 
in  the  literature,  much  less  a  theoretical  prediction  of  such  a 
relationship . 

The  following  chapters  contain  discussions  of  the  experimental 
approach  and  observed  results .  These  results  will  be  interpreted  in 
terms  of  the  conclusions  drawn  from  Equations  (2.19)  and  (2.27). 


CHAPTER  III 


EXPERIMENTAL  PROCEDURES 

At  the  time  that  the  present  investigation  was  begun,  very  little 
research  held  been  performed  and  reported  on  NAI.  A  program  was 
therefore  initiated  and  executed  with  the  purpose  of  producing  NAI 
in  certain  propellants  that  were  readily  available  at  the  University 
of  Utah.  The  initial  experiments  made  use  of  shaped  grains  whose 
burning  area  changed  (accompanied  by  a  corresponding  change  in  the 
chamber  pressure)  as  burning  progressed.  It  was  hoped  that  this 
pressure  perturbation  would  be  adequate  to  trigger  instability. 

Although  many  variations  were  employed  by  using  variously  shaped  grains, 
neither  oscillatory  instability  nor  chuffing  was  observed  in  some  fifty 
runs . 

It  was  later  discovered  that  a  few  of  these  runs  were  actually 
made  in  what  we  now  know  is  the  unstable  regime.  The  fact  that  instability 
was  not  detected  was  probably  due  to  the  manner  in  which  the  pressure 
transducer  was  mounted.  The  transducer  was  mounted  inside  a  Kistler, 

WC-l4,  water-cooled  adaptor  located  in  the  wall  of  the  burner.  The 
transducer  vas  separated  from  the  gas  chamber  of  the  burner  by  a 
column  of  silicone  grease  approximately  l/8  inch  in  diameter  and  3/4 
inch  long  for  the  runs  in  question.  When  this  adaptor  was  utilized 
in  connection  with  the  low-L*  burner  (in  which  NAI  is  readily  observed) 
oscillations  were  apparently  damped  out  by  the  column  of  grease. 
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Therefore,  the  results  of  the  experiments  with  the  shaped  grains  are 
considered  a?  Inconclusive.  However,  these  results  did  lead  to  the 
decision  to  design  a  burner  that  would  operate  at  very  short  stay  times 
(i.e.,  low  L*). 

A.  The  Low-L*  Burner  and  Its  Instrumentation 

The  low-L*  burner  that  was  designed  was  modeled  after  one  used 
at  the  Naval  Ordnance  Test  Station  [26,28].  Figure  4  shows  a  schematic 
drawing  of  the  assembled  burner.  Hie  burner  consists  of  two  10. 5- 
inch  flanges  with  a  short  section  of  1.5-lnch  i.d.,  heavy  walled  pipe 
between  the  f  langes .  In  order  to  be  able  to  vary  the  length  of  the 
burner,  several  sections  of  the  pipe  ranging  from  one-half  Inch  to 
3  inches  In  length  were  made. 

A  cylindrical  disc  of  propellant  1.5  inches  In  diameter  and 
usually  one  quarter  to  three  quarters  inch  thick  was  placed  In  the 
pipe.  This  section  was  then  positioned  between  the  two  flanges  so 
that  the  propellant  was  flush  against  a  flange  containing  a  Kistler 
Model  hOl  pressure  transducer.  A  carbon  nozzle  1  Inch  In  diameter 
controlled  the  flow  of  exhaust  gases  and  was  situated  in  the  flange 
opposite  the  propellant.  Each  nozzle  was  machined  so  that  it  would 
fit  flush  with  the  surface  of  the  flange  as  shown  in  Figure  k.  Figure 
5  is  a  photograph  of  three  of  the  actual  nozzles  used  in  the  experiments. 
The  nozzles  were  made  with  throat  diameters  ranging  from  0.180  inches 
to  0.328  inches.  This  resulted  in  a  range  of  Kn  (the  ratio  of  burning 
area  to  nozzle  throat  area)  between  21  and  69  and  allowed  L*  to  be 
varied  (for  burner  lengths  of  3  inches  or  less)  from  approximately 
5  cm.  to  well  over  300  cm. 
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LOW  L*  MOTOR 


Figure  4.  Schematic  drawing  of  lov-L*  burner. 
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Figure  5.  Graphite  flow  control  nozzles  used  in  the  low-L* 
burner. 
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The  pressure  signal  from  the  transducer  passed  through  a  Kistler 
Model  568  charge  amplifier,  was  displayed  on  a  Type  535  Tektronix 
oscilloscope  where  it  was  photographed  using  Polaroid  film,  and  was 
recorded  on  a  Model  207,  Precision  Instrument  magnetic  tape  recorder. 

For  convenience,  the  signal  was  also  observed  visually  on  a  Type  150A 
Hewlett-Packard  oscilloscope.  The  tape  recorder  was  operated  at  a 
tape  speed  of  15  inches  per  second  and  played  back  at  7  and  one-half 
inches  per  second  giving  a  two-to-one  expansion  of  the  time  Beale. 

B.  The  Low-L*  Burner  Adapted  for  High  Speed  Motion  Pictures 

During  the  course  of  the  experimental  studies  it  was  felt  that 
high-speed  motion  pictures  of  the  unstable  combustion  process  might  be 
informative.  In  order  to  facilitate  this,  plexiglas  windows  were  con¬ 
structed  with  flow  control  nozzles  machined  at  their  center.  The 
plexiglas  was  then  attached  to  the  low-L*  burner  in  the  manner  shown 
in  Figure  6 .  A  Wollensak,  Type  WF3T,  Fastax  high  speed  camera  using 
sometimes  black  and  white  and  sometimes  color  film  (Eastman  Kodak  Plus-X 
and  Ektac hrome - ER  film)  was  operated  at  speeds  between  500  and  1000 
frames  per  second  in  making  the  movies . 

The  results  from  these  studies  were  interesting  but  not  quantitative. 
They  showed  that  the  light  intensity  from  the  burning  propellant  oscil¬ 
lated  in  the  same  way  that  the  pressure  does,  but  due  to  the  fact  that 
the  camera  and  the  oscilloscope  were  operated  independently  of  each 
other,  a  quantitative  measurement  of  this  effect  could  not  be  made. 


(ADAPTED  FOR  MOVIES) 


Figure  6.  Low-L*  burner  adapted  for  movies. 
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C.  The  Vacuum  System 

After  many  experimental  runs  had  been  made,  it  was  observed  that 
MAI  appears  to  occur  more  readily  at  low  operating  pressures  [see  also 
2,20,28].  A  disadvantage  of  operating  at  low  pressures  is  that  if  the 
ratio  of  the  chamber  pressure  to  the  ambient  press  ore  falls  below  a 
certain  value,  critical  flow  will  not  be  maintained  in  the  nozzle.  The 
exact  relationship  may  be  derived  as  follows:  the  stagnation  tempera¬ 
ture  for  sonic,  isentropic  flow  is 

T  =  T  (3.1) 

ch  t  v  2  ' 

where  T  is  the  ratio  of  specific  heats.  Combining  Equation  ( 3 •  1 ) 
with  the  isentropic  relationship  between  pressure  and  temperature 
yields 

r 

fch  „  (i_±-£)  r  -1  (3-2) 

For  the  propellants  used  in  this  study  the  ratio  of  specific  heats 
ranges  from  1.20  to  1.24  (see  Appendix  B).  Substituting  these  values 
of  r  into  (3.2)  yields  a  critical  pressure  ratio  of  approximately  1.8. 
This  means  that  if  a  ratio  of  approximately  two  is  not  maintained 
between  the  chamber  pressure  and  the  ambient  pressure,  then  critical 
flow  does  not  exist  in  the  nozzle. 

Many  r\ins  were  made  with  sub-critical  pressure  ratios.  The 
results  are  therefore  in  doubt,  because  their  interpretation  requires 
critical  nozzle  flow.  In  order  to  assure  critical  flow  for  low  pres¬ 
tank  was  constructed  into  which  the  experimental 


sure  firings  a  vacuum 
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burner  could  exhaust  at  sub-atmospheric  pressures.  Use  of  the  vacuum 
+.pwv  for  this  purpose  made  possible  the  investigation  of  chamber  pres¬ 
sures  slightly  above  atmospheric.  Although  sub -atmospheric  pressures 
could  also  have  been  investigated,  it  was  felt  that  this  would  have 
been  beyond  the  scope  of  the  present  study.  The  vacuum  tank  was  con¬ 
structed  by  the  Lang  Wayne  Equipment  Company  of  Salt  Lake  City,  Utah. 

It  45  cubic  foot  volume  and  will  withstand  a  vacuum.  A  Type  MD 
674,  Nash  Hytor  Vacuum  Pump  driven  by  a  5-horsepower  motor  was  utilized 
in  series  with  a  No.  2-26-6  Nash  Air  Ejector  to  draw  the  necessary 
vacuum  on  the  tank.  Figure  7  shows  a  schematic  drawing  of  the  vacuum 
system. 

The  tank  was  constructed  so  that  the  experimental  burner  could  be 
attached  directly  to  it.  A  4-inch  i.d.  pipe,  12  inches  long  screwed 
into  the  access  opening  of  the  tank,  and  a  10.5-inch  flange  was  welded 
onto  the  free  end  of  the  pipe  so  that  the  burner  could  be  bolted  directly 
to  the  flange.  Two  wires  were  cemented  into  drilled  holes  in  the  wall 
of  the  access  pipe  to  provide  pres sure -tight  ignition  leads.  Figure  8 
is  a  photograph  of  the  burner  attached  to  the  vacuum  tank.  The  ignition 
leads  can  be  seen  at  the  top  of  the  photograph,  and  the  charge  amplifier 
can  be  seen  alongside  the  burner.  The  tank  itself  was  placed  outside 
the  building  and  the  access  pipe  passed  through  a  wall  in  such  a  manner 
that  the  flange  was  indoors . 

The  effect  of  exhausting  to  a  pressure  lower  than  atmospheric 
can  best  be  demonstrated  by  the  comparison  of  the  pressure  traces 
of  rums  made  under  similar  conditions,  with  and  without  the  vacuum 


1  .  VACUUM  TANK 

2.  ASSEMBLED  LOW-lT  BURNER 

3.  BLEED  VALVE 

4.  ABSOLUTE  PRESSURE  GAUGE 

5.  NH3  TANK 

6.  NASH  NO.  2-26-6  AIR  EJECTOR 

7.  NASH  HYTOR  VACUUM  PUMP,  TYPE  MD674 


Figure  7-  Schematic  diagram  of  the  vacuum  system. 
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Figure  8.  Assembled  low-L*  burner  attached  to  vacuum 
tank.  Ignition  wires  are  connected  and  the  charge  amplifier 
can  also  be  seen. 


tank.  Figure  9  shows  a  comparison  of  two  such  runs.  Although  the 
time  scale  is  the  same  for  the  two  traces,  the  bottom  picture  does 
not  include  the  total  run  due  to  the  slow,  irregular  nature  of  the 
chuffing.  The  general  effect  of  exhausting  into  a  vacuum  was,  as  can 
be  seen:  delay  times  of  several  seconds  often  occurring  between  chuffs 
with  a  few  oscillations  taking  place  on  the  chuffs  before  extinction. 
However,  when  exhausting  to  atmospheric  pressure,  rapid-fire  chuffing 
occurs  with  few,  if  any,  oscillations  being  developed  on  a  chuff .  Another 
effect  of  utilizing  the  vacuum  tank  was  that  the  propellant  would  often 
extinguish  without  reignition  when  the  ambient  pressure  was  less  than 
5  psia,  whereas  at  atmospheric  pressure  samples  were  seldom  completely 
extinguished. 

D.  Propellants 

The  initial  experiments  that  were  carried  out  with  shaped  gralfas, 
utilized  G  and  F  propellants.  The  compositions  and  mixing  procedure 
for  these  propellants  were  established  at  the  University  of  Utah  by 
Ralph  Coates  [8]  and  are  recorded  in  Appendix  B  along  with  their 
physical  properties.  When  positive  results  were  not  obtained  from 
these  initial  experiments,  it  was  suggested  [27]  that  a  new  burner  be 
designed  and  that  a  polyurethane  propellant  be  used.  As  a  result, 
the  burner  described  in  Section  A  of  this  chapter  was  designed  and 
constructed.  A  polyurethane  propellant  was  made  and  burned  in  the  low-L* 
burner  giving  positive  results.  The  formulation  which  was  used  was 
ob-ained  from  the  Naval  Ordnance  Test  Station  and  is  recorded  in 
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Run  No.  310-233,  L*  -  30  cm.,  L*  =  56  cm. 
Ambient  pressure:  atmospheric 


Run  No.  49-181.  L*  =  38  cm.,  L*  =  56  cm. 
Ambient  pressure:  5  psia 

(This  is  a  play-back  from  the  tape  recorder.) 


Figure  9*  A  comparison  of  pressure- time  traces  for  similar  runs 
made  with  and  without  the  vacuum  tank.  Both  runs  were  made  with  F 
propellant.  The  final  burn-out  is  not  shown  in  the  bottom  trace. 

The  sweep  rate  is  1  second/major  division,  and  the  pressure  scale  is 
10  psi/major  division  for  both  traces. 
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The  extremely  high  viscosity  of  the  propellant  in  the  uncured 
state  discouraged  extensive  use  of  this  propellant  due  to  the  prob¬ 
lems  encountered  in  the  mixing  and  casting  process.  Therefore,  after 
approximately  forty  runs  had  been  made,  it  vas  decided  to  discontinue 
use  of  the  propellant  even  though  positive  results  had  been  obtained. 
The  raw  data  from  these  runs  can  bn  found  in  Appendix  E,  Table  4. 

The  experience  gained  from  the  runs  made  with  the  polyurethane  pro¬ 
pellant  provided  the  experimental  understanding  necessary  to  obtain 
positive  results  with  G  propellant  by  operating  at  very  low  pressures 
and  low  values  of  L*.  Subsequently,  NAI  was  also  observed  in  F  pro¬ 
pellant  even  though  other  investigators  had  not  been  able  to  observe 
NAI  in  comparable  propellants . 

From  the  above  observations  it  was  decided  to  use  the  PBAA 
propellants  (F  and  G)  in  the  present  investigation.  A  third  propel¬ 
lant  containing  five  per  cent  aluminum  (displacing  perchlorate  from 
F  propellant)  was  formulated  and  designated  TF  propellant.  These 
three  propellants  formed  the  basis  of  the  investigation.  Exploratory 
work  was  initiated  with  two  other  propellants  designated  GB  and  XF. 

The  GB  propellant  is  an  uncatalyzed  propellant  similar  to  G  excepting 
that  it  contains  two  per  cent  carbon  black.  The  purpose  of  the  carbon 
black  is  to  reduce  the  amount  of  radiation  penetrating  the  solid 
(G  propellant  is  slightly  translucent).  The  XF  propellant  contains 
ten  per  cent  aluminum  and  is,  therefore,  more  highly  aluminized  than 
any  of  the  other  propellants.  Appendix  B  contains  the  formulations 
and  physical  properties  of  all  the  propellants  with  the  exception  of 


the  physical  properties  of  the  polyurethane  propellant.  Virtually  no 
conclusions  were  drawn  from  the  use  of  this  propellant  and,  therefore, 
its  properties  have  not  been  recorded. 

E.  Experimental  Run  Procedure 

In  the  initial  experiments  with  shaped  grains  a  great  deal  of 
time  was  spent  in  preparing  the  samples  for  firing.  The  samples  were 
machined  to  the  proper  size,  and  then  several  coats  of  a  burning 
restricter  were  applied  to  the  side6  of  the  sample  to  insure  even 
burning.  The  second  experimental  apparatus  (the  low-L*  burner)  was 
constructed  taking  into  consideration  all  of  the  disadvantages  of  the 
first.  The  design  of  the  burner  is  very  simple,  and  this  has  resulted 
in  a  minimum  amount  of  effort  in  preparing  a  sample  for  firing.  It 
was  hoped  that  the  simplicity  Of  the  experimental  system  might  help 
to  compensate  for  the  complexity  of  the  physical  process. 

The  propellant  was  cast  into  the  same  heavy  wall  tubing  that  was 
used  in  the  construction  of  the  walls  of  the  burner.  However,  the 
molds  were  reamed  out  slightly  in  order  to  allow  for  shrinkage  as  the 
propellant  cooled.  The  samples  fit  snugly  in  the  burner.  In  order  to 
prepare  a  sample  for  a  firing,  a  piece  of  propellant  approximately  the 
desired  length  was  cut  from  the  cast  grain  with  a  hacksaw.  The  grains 
were  usually  fifteen  to  eighteen  inches  long.  The  surface  of  the 
sample  was  then  scraped  smooth  with  a  razor  blade  until  it  formed  a 
right  angle  with  the  side  of  the  sample  and  was  the  desired  length. 

The  sample  was  then  forced  into  the  test  section  that  was  to  be  used. 
In  some  of  the  experiments  the  test  section  was  coated  lightly  with  a 
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a  silicone  grease  in  order  to  facilitate  the  entry  of  the  sample 
into  the  test  section.  Later,  as  an  extra  precaution  against  movement 
of  the  sample  during  the  firing, or  "burning  down  the  side,  the  sample 
was  cemented  to  the  walls  of  the  test  section  with  an  epoxy  resin. 

No  appreciable  difference  was  noted  between  runs  made  using  the  dif¬ 
ferent  procedures. 

Samples  that  did  not  ignite  uniformly  would  burn  irregularly 
causing  an  uneven  pressure  trace,  G  propellant  being  the  most  difficult 
of  the  propellants  to  ignite.  It  was  also  the  slowest  burning  pro¬ 
pellant  and  was,  therefore,  more  susceptible  to  uneven  burning.  It 
should  be  noted  that  the  requirement  for  the  strength  of  the  ignition 
pills e  also  varied  with  the  length  of  the  test  section.  For  these 
reasons  the  ignition  process  had  to  be  tailored  to  fit  the  needs  of 
almost  each  separate  run,  particularly  when  the  sample  to  be  burned 
was  of  G  propellant. 

The  primary  source  of  ignition  consisted  of  a  length  of  an 
exploding,  exothermic  fuse  wire  (Lo-R  Pyrofuze  made  by  the  Pyrofuze 
Corp.).  The  wire  explodes  giving  off  heat  when  its  temperature 
increases  sufficiently  due  to  the  heating  from  a  current  passing 
through  the  wire.  A  voltage  drop  of  approximately  60  to  70  volts  was 
necessary  to  ignite  the  Pyrofuze.  For  a  typical  run  involving  G 
propellant  a  „vo-inch  length  of  0.015-inch  diameter  Pyrofuze  wire  was 
used.  With  the  other  propellants  a  length  approximately  1.5  inches 
long  of  0.010-inch  diameter  wire  was  used.  A  pyrotechnic  paste  com¬ 
posed  of  finely  powdered  boron,  aluminum,  and  ground  crystals  of 
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potassium  perchlorate  with  polyisobutylene  as  a  carrier,  was  painted 
on  tne  surface  of  the  sample  as  a  secondary  source  of  ignition  when 
test  conditions  warranted  it.  The  paste  was  dissolved  in  methylene 
chloride  which  evaporated  when  the  paste  was  spread  on  the  sample. 

After  the  propellant  sample  and  the  ignitor  had  been  prepared, 
the  following  procedure  was  carried  out:  a  thin  coating  of  silicone 
grease  was  spread  over  the  transducer  to  protect  it  as  the  propellant 
burned  out,  and  the  test  section  was  placed  between  the  flanges 
with  the  ignitor  lead  wires  extending  from  the  nozzle.  The  flanges 
were  then  securely  bolted  together.  For  the  runs  that  utilized  the 
vacuum  tank,  the  necessary  vacuum  was  drawn  and  the  system  sealed  off 
(the  vacuum  pump  was  not  operated  during  a  firing).  The  tank  was 
not  completely  evacuated  for  each  run  because  the  ignition  problems 
became  more  acute  as  the  pressure  level  became  lower.  The  pressure 
was  maintained  only  at  a  level  sufficiently  low  to  ensure  critical  flow 
in  the  nozzle  during  the  run.  The  sample  was  then  ignited  and  the  pres¬ 
sure  trace  recorded .  After  the  completion  of  a  run,  a  small  amount  of 
ammonia  was  bled  into  the  vacuum  tank  in  order  to  neutralize  the  acidic 
combustion  gases  and  their  deposits .  The  tank  was  then  flushed  out  well 
with  compressed  air.  After  several  runs  the  tank  was  washed  down  with 
water  and  allowed,  to  dry  before  continued  use. 


CHAPTER  IV 


ANALYSIS  OF  EXPERIMENTAL  RESULTS 

The  results  of  the  experimental  firings  that  were  made  in  con¬ 
nection  with  the  present  study  have  been  recorded  on  magnetic  tape 
as  pi-es  sure -time  traces.  The  present  chapter  discusses  the  measure¬ 
ments  made  from  these  traces  and  the  conditions  under  which  the 
firings  were  made.  The  tabulation  of  these  resultB  has  been  divided 
into  three  sets  of  tables.  The  first  set  of  tables  contain  what 
can  be  termed  "raw  data"  (i.e.,  the  pertinent  physical  measurements). 
The  second  and  third  sets  include  data  that  has  been  analyzed  by  two 
different  methods.  The  first  part  of  this  chapter  discusses  these 
tables,  and  the  remainder  of  the  chapter  is  devoted  to  a  discussion 
of  the  criteria  used  in  analyzing  the  data,  concluding  with  a  series 
of  graphs  containing  the  data. 

A.  A  General  Discussion  of  the  Reported  Data 

Tables  4-9  contain  the  raw  data  pertinent  to  each  run,  including 
the  nozzle  throat  diametei,  the  initial  length  of  the  propellant 
samples  (Lp),  the  chamber  length  (L,h),  and  the  ambient  pressure. 
Observations  of  the  general  behavior  of  the  runs  are  included  as  the 
range  of  the  mean  pressure  (the  mean  pressure  during  a  firing  was 
not  always  constant  particularly  during  unstable  combustion),  and  a 
brief  note  as  to  whether  the  firing  was  stable  or  unstable.  The 
code  which  was  used  for  the  latter  entry  Is  quite  simple:  C  stands 
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for  chuffing;  0  stands  for  oscillatory  burning;  S  stands  for  stable 
burning;  £  means  that  the  propellant  extinguished;  and  U  means  that 
the  pressure  vas  uneven  or  tuat  a  form  of  irregular  burning  was 
apparently  present.  Each  separate  table  contains  data  for  one  pro¬ 
pellant  only,  and  the  Individual  runs  are  grouped  according  to  nozzle 
size  beginning  with  the  lowest  value  of  L*  and  increasing. 

The  second  set  of  tables  contains  information  concerning  the 
stability  limit  as  a  function  of  L*  and  pressure.  This  limit  can  be 
defined  mathematically  as  the  set  of  conditions  at  which  the  growth 
factor  (tj)  equals  zero  or  in  experimental  terms,  as  the  point  where 
the  pressure  oscillations  begin  to  decay  rather  than  grow.  Therefore, 
during  the  course  of  a  firing,  if  chuffing  or  oscillatory  burning  is 
followed  by  stable  burning,  then  the  stability  limit  has  been  crossed, 
and  the  value  of  L*  at  this  limit  is  the  critical  L*.  Tables  lO-?  3 
contain  the  tabulated  values  that  were  determined  to  be  the  critical 
L*  for  the  various  propellants  tested  and  the  pressure  at  which  the 
evaluation  was  made.  Section  D  of  this  chapter  describes  the  method 
used  to  determine  values  of  L*  from  the  pressure-time  traces.  The 
data  from  these  tables  have  been  plotted  (see  Figure  10)  as  suggested 
by  the  analysis  of  Akiba  and  Tanno  [l],  and  Sehgal  and  Strand  [35] 
showing  the  stability  limits  for  the  four  propellants,  TF,  XF,  F,  and 
0.  In  order  to  eliminate  the  confusion  that  would  result  from  plotting 
all  of  the  data  points  of  the  different  propellants,  only  the  lines 
representing  the  data  are  shown  in  the  figure.  The  figure  indicates 
that  the  aluminized  propellants,  XF  and  TF,  are  more  unstable  than 
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Figure  10.  Stability  limits  of  the  propellants  tested  in  this 
study.  Plots  represent  the  best  line  drawn  through  the  data  of 
Tables  10-13,  Appendix  F. 


-  51  - 


the  non- aluminized  propellants,  and  that  the  slower  burning,  uncatalyzed 
G  propellant  is  more  stable  than  the  catalyzed  propellants. 

Tables  14-18  contain  the  tabulated  results  of  the  frequency 
measurements,  the  observed  pressure,  and  the  L*  calculated  at  the 
point  where  the  frequency  was  measured.  The  dimensionless  frequency, 

7,  and  the  dimensionless  L*,  £,  are  also  tabulated  in  these  tables. 

L*- frequency  data  were  not  obtained  from  many  firings  which  exhibited 
chuffing  or  oscillatory  burning  because  they  were  not  recorded  on  the 
tape  recorder,  or  because  the  recorder  did  not  function  properly. 

Data  from  traces  that  showed  evidence  of  uneven  burning  or  whose 
pressure  varied  greatly  were  not  included  in  these  tables  even  when 
they  were  available  because  of  the  irregularities  involved. 

The  steady-state  burning  rate  enters  into  the  expressions  for 
both  the  dimensionless  frequency  and  L*.  The  value  that  was  used  in 
calculating  these  dimensionless  variables  was  determined  in  the  following 
manner:  the  mean  pressure  of  the  oscillations  was  measured  as  accurately 
as  possible  from  the  pressure  trace,  and  the  linear  burning  rate  was 
then  evaluated  for  the  measured  pressure  from  data  that  had  been  ob¬ 
tained  previously  in  strand  bomb  tests.  The  burning  rate  curves  for 
the  propellants  tested  are  included  in  Appendix  B  with  the  physical 
properties . 

Table  14  also  contains  a  tabulation  of  the  growth  factor,  tj, 
that  was  determined  for  some  of  the  firings  involving  TF  propellant. 

The  manner  of  approximating  the  values  of  this  variable  will  be  explained 
in  Section  D  of  this  chapter. 
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B.  Discussion  of  Typical  Pressure-Time  Traces 

The  photographs  of  two  pressure-time  traces  that  show  the  dif¬ 
ferent  aspects  of  non-acoustic  instability  quite  clearly  make  up 
Figure  11.  The  two  traces  are  from  firings  made  vith  7  and  TF  pro¬ 
pellant  burning  under  approximately  the  same  conditions.  The  first 
chuff  in  both  traces  was  of  short  duration  and  the  second  was  con¬ 
siderably  longer  with  an  oscillating  pressure  on  the  crest  of  the 
chuff.  This  indicates  that  the  conditions  are  in  a  slightly  more 
stable  regime  which  is  in  agreement  with  the  observation  that  stability 
improves  with  increasing  L*  (L*  increases  as  burning  progresses).  After 
chuffing  twice,  both  samples  reignited  and  remained  burning  even  though 
they  passed  through  a  period  of  severe  oscillations.  From  the  previous 
discussion  it  can  be  understood  that  the  point  at  which  the  amplitude 
of  the  oscillating  pressure  is  a  maximum  defines  the  stability  limit. 

The  oscillations  that  continued  beyond  this  point  were  decreasing  in 
amplitude  and  could,  therefore,  be  considered  as  occurring  in  a  stable 
region.  This  implies  that  a  small  perturbation  in  the  pressure  would 
not  grow  but  would  be  suppressed,  and,  presumably,  that  a  large  dis¬ 
turbance  would  decay  also.  Toward  the  end  of  the  bottom  trace  such 
a  disturbance  occurred  causing  an  immediate  increase  in  the  oscillating 
pressure  (i.e.,  the  amplitude  did  not  grow  slowly  as  at  the  first  of 
the  trace).  This  phenomenon  has  been  interpreted  as  being  due  to  an 
air  bubble  in  the  propellant  or  something  Bimilar  that  would  cause  an 
almost  instantaneous  pressure  disturbance. 

The  rapid-fire  pressure  oscillations  at  the  end  of  the  trace 
are  attributed  to  uneven  residual  burning.  The  reduced  burning  area 


Run  No.  410-163 
L*  range:  8.6  -  46  cm. 
Ambient  pressure:  20  psla 


F  Propellant 
P  ■  70  psla 


Run  No.  48-011  TF  Propellant 

L*  range:  2.0  -  40.6  cm.  P  »  78  psla 
Ambient  pressure:  12.5  psla 

Figure  11.  Typical  pressure-time  traces  demonstrating  NAI . 
The  ordinate  for  both  photographs  is  20  psl/major  division  and  the 
time  scale  is  0.25  sec. /major  division. 


causes  the  pressure  to  decrease  until  it  is  low  enough  to  enter  an 
unstable  regime  for  the  value  of  L*  at  those  conditions. 

It  should  be  pointed  out  here  that  the  frequency  of  the  oscil¬ 
lations  varied  throughout  the  period  of  instability  in  both  traces . 

This  variation  in  frequency  for  changing  L*  at  constant  pressure  was  a 
general  observation  for  the  'ings  made  with  F,  IF,  and  XF  propellants. 

C .  The  Anomalous  Behavior  of  G  Propellant 

The  photograph  of  the  trace  in  Figure  12  corresponds  to  a  firing 
made  with  G  propellant,  and  is  quite  similar  in  its  general  nature  to 
those  in  Figure  11.  A  difference  between  this  trace  and  the  two 
previous  is  that  the  oscillations  attain  a  "pseudo  steady  state"  where 
they  maintain  a  constant  amplitude,  neither  growing  further  nor  decaying. 
This  phenomenon  is  typical  of  firings  made  with  G  propellant.  The 
oscillations  would  reach  a  certain  amplitude  and  then  continue  at 
that  amplitude  for  a  considerable  time.  Apparently,  losses  for  the 
slower  burning,  uncatalyzed  propellant  were  great  enough  to  inhibit 
further  growth  in  the  amplitude  of  the  pressure  oscillations.  This 
type  of  behavior  was  not  observed  in  firings  made  with  any  of  the 
catalyzed  propellants . 

Another  significant  difference  between  this  trace  and  the  other 
two  is  that  the  frequency  changed  very  little  during  the  period  of 
constant  amplitude  oscillations,  whereas  with  the  other  propellants, 
the  frequency  changed  continually  during  a  run.  This  observation  was 
also  found  to  be  general  for  G  propellant,  whereas  no  comparable 
behavior  was  observed  with  any  of  the  other  propellants.  Another 


Run  Number 
L*  range: 
Top  trace: 

Bottom  trace: 
Sweep  rate: 


39-233  G  Propellant 

13  to  8l  cm.  P  =  23  to  30  psla 
Oscillating  component  of  pressure, 

5  psi/major  division. 

Total  pressure  10  psi/major  division 
2  seconds/major  division 


Figure  12.  Pres sure- time  trace  for  G  propellant.  Note  the 
constant  amplitude  oscillations  and  the  slight  change  in  this  amplitude 
near  the  first  part  of  the  firing. 
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observation  that  can  be  made  of  the  trace  in  Figure  12  concerns  the 
slight  change  in  amplitude  that  is  observed  shortly  after  the  initiation 
of  the  constant  amplitude  oscillations .  Examination  of  the  expanded 
trace  reveals  the  presence  of  vhat  appears  to  be  "beating"  or  the 
existence  of  oscillations  of  different  frequencies.  Ibis  phenomenon 
was  observed  in  three  or  four  runs  where  the  physical  conditions  were 
essentially  identical. 

From  the  above  comparison  of  the  firings  made  with  the  three 
propellants  (F,  TF,  and  G),  it  appears  reasonable  to  postulate  that 
the  mechanism  causing  instability  i n  the  same  for  TF  and  F  propellants 
but  is  probably  different  for  G  propellant.  An  investigation  of  the 
properties  of  the  propellants  in  an  attempt  to  explain  the  reason  for 
this  difference  in  behavior  shows  that  the  burning  rate  of  G  propellant 
is  approximately  half  that  of  the  other  propellants  at  comparable 
pressures  (see  Figure  29,  Appendix  B).  A  second  difference  in  the  prop¬ 
erties  of  the  propellents  is  involved  in  their  translucence.  G  propel¬ 
lant  is  a  very  light  tan  color  and  is  slightly  translucent  to  visible 
light,  while  TF  and  F  propellants  both  contain  the  black-colored  copper 
chromite  causing  them  to  be  opaque. 

In  an  attempt  to  determine  which  of  these  differences  had  the 
greatest  effect  on  the  behavior  of  G  propellant,  another  propellant, 
designated  GB,  was  compounded  having  approximately  the  same  burning 
rate  as  G  but  that  was  opaque  to  visible  light.  The  opaqueness 
resulted  from  adding  two  per  cent  carbon  black  (replacing  ammonium 
perchlorate)  to  the  composition  of  G  propellant.  The  results  of  the 
firings  made  with  this  propellant  comprise  Tables  9  and  18,  and  show 
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that  it  did  behave  differently  from  G  propellant  when  fired  in  similar 
conditions.  A  more  complete  discussion  of  the  results  obtained  from 
these  two  propellants  will  be  found  later  in  this  chapter. 

D.  Analysis  of  Data 

The  quantitative  data  that  were  gleaned  from  the  pressure  traces 
included  the  frequency,  the  L*  at  locations  where  the  frequency  wes 
measured,  the  critical  L*  where  applicable,  the  pressure  at  those  times 
when  other  measurements  were  made,  and  for  certain  runs  involving  TF 
propellant,  the  growth  factor  .  Although  these  quantities  are  not 
specifically  raw  data,  they  were  measured  or  calculated  directly  from 
the  pressure  traces  and  will  be  referred  to  as  data  throughout  the  remainder 
of  the  text.  The  frequency  was  measured  from  the  expanded  traces  that 
were  obtained  from  the  play-back  of  the  permanent  record  made  with  the 
tape  recorder.  Figure  3  is  a  typical  example  of  the  type  of  record 

that  was  used  for  measuring  the  frequency.  Many  of  the  nans  contained 

chuffs  with  only  four  or  five  cycles  of  oscillations  whose  amplitude 
was  varying  widely.  For  data  of  this  nature,  the  frequency  was  measured 
over  three,  four,  and  five  cycles  and  then  averaged. 

In  order  to  correlate  frequency  and  L*  it  was  necessary  to  deter¬ 
mine  the  value  of  L*  at  the  same  location  on  the  trace  where  the  frequency 

was  measured.  This  was  done  by  using  a  planimeter  to  integrate  the 
pressure-time  trace.  The  area  beneath  the  curve  of  an  oscilloscope 
trace  approximately  represents  the  mass  of  propellant  burned,  and  the 
change  in  L*  is  directly  proportional  to  the  mews  of  propellant  burned, 
ewsuming  that  the  burning  area  is  constant  throughout  a  firing.  This 
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method  of  determining  values  of  L*  evolves  from  considering  equilibrium 
mass  flow  through  a  sonic  nozzle.  The  resulting  equation  is 

m  -  j Pdt  (^-1) 

or, 

dt  (^-2) 

which  is  to  say  that  the  increase  in  L*  over  the  initial  L*  is  pro¬ 
portional  to  the  area  beneath  the  pressure-time  trace.  The  values  of 
L*  at  the  locations  in  question  were  determined  by  integrating  the 
traces  with  a  planimeter  and  averaging  three  such  readings. 

It  has  been  mentioned  previously  that  the  growth  factor  was 
approximated  for  many  of  the  firings  involving  TF  propellant.  A 
standard  method  of  determining  coefficients  in  exponential  terms  was 
utilized  in  calculating  the  values  of  q.  The  double  anqplitude  of  the 
pressure  oscillations  was  plotted  verus  time  on  semi-log  coordinates. 
The  measured  slope  of  the  data  plotted  in  this  manner  is  the  growth 
factor. 

A  simple,  theoretical  method  of  correlating  the  growth  factor 
with  the  frequency  and  L*  is  not  available  at  the  present.  Therefore, 
a  trial- and- error  method  was  employed  to  determine  an  empirical 
relationship  between  these  variables .  From  this  study  it  appeared 
that  the  best  correlation  of  the  data  was  obtained  by  plotting  the 
growth  factor  versus  the  reciprocal  of  the  product  of  the  dimensionless 
frequency  and  the  dimensionless  L*  to  the  one -half  power.  Figure  13 
represents  data  for  TF  propellant  plotted  in  this  manner.  Although  the 
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best  line  through  the  data  does  appear  to  be  curved,  for  simplicity  a 
straight  line  through  76^^  equal  to  two  could  also  be  considered  a 
reasonable  representation  of  the  data.  The  results  obtained  from  this 
figure  and  their  application  will  be  discussed  more  in  detail  in  the 
next  chapter. 

The  majority  of  the  growth  factor  values  were  obtained  from 
series  of  ten  or  fewer  pressure  oscillations  with  smooth  curves  drawn 
through  the  pressure  peaks  to  define  envelopes.  (See  Figure  3  for  an 
example  of  the  type  of  expanded  pressure  traces  used.)  The  precision 
leaves  much  to  be  desired,  but  is  probably  better  than  implied  by  the 
heuristic  correlation  of  Fig.  13 • 

E.  Presentation  of  Experimental  Results 

There  are  several  methods  by  which  the  data  obtained  in  this 
study  can  be  correlated.  Figure  10,  containing  the  data  analyzed  in 
terms  of  a  stability  boundary,  bas  already  been  discussed.  Data  have 
also  been  plotted  according  to  the  preferred  f requency-pressure 
relationship  that  was  discussed  briefly  in  Chapter  I.  Figure  14  con¬ 
tains  four  plots  for  the  five  propellants,  F,  TF,  XF,  0,  and  GB.  As 
predicted  by  the  preferred  frequency  concept,  the  data  fall  into  fairly 
narrow  band  of  frequencies  for  a  given  pressure.  The  data  obtained 
with  the  catalyzed  propellants  appear  to  be  consistent  and  not  scattered 
greatly,  but  the  data  from  the  uncatalyzed  propellants  (particularly  G) 
do  not  show  the  same  type  of  trend  and  are  quite  scattered.  For  compara¬ 
tive  purposes,  data  reported  by  Eisel,  et  al.  [12]  for  a  highly  aluminized 
PBAA  propellant  have  been  included  with  the  data  for  XF  propellant.  The 
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Figure  14.  Preferred-frequency  pressure  plot  of  data. 
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from  the  two  propellants  are  In  excellent  agreement  even  though 
aluminum  content  is  probably  higher  for  Eisel's  propellant  than 
XF  propellant.  It  appears  that  this  method  of  correlating  instability 
.a  has  merit. 

It  has  been  noted  previously  that  the  frequency  of  the  oscillations 
xied  during  a  firing  as  L*  varied,  and  the  values  of  L*  and  the  fre- 
aencies  have  been  determined  for  these  conditions.  Data  of  this 
ature  appeared  to  correlate  well  when  plotted  as  frequency  versus  the 
•eciprocal  of  L*,  and  those  data  obtained  for  TF  propellant  have 
therefore  been  presented  as  Figure  15*  The  various  nozzle  sizes, 
corresponding  to  different  symbols  on  the  curve,  represent  different 
pressures,  and  a  pressure  effect  is  apparent.  This  effect  was  elimi¬ 
nated  by  using  the  dimensionless  variables  as  suggested  by  the  form  of 
Equation  (2.19).  This  equation  can  be  written  as 


7 


1  F  tan  £  -  1 _ 

!  _F  (1-T}tan£)  +  T]  +  tan  £ 


(2.19) 


The  form  of  this  equation  suggests  that  7  be  plotted  versus  l/£.  Figures 
16,  17,  and  18  contain  the  data  of  TF,  F,  and  XF  propellants,  respectively, 
plotted  in  this  manner.  The  same  symbols  have  been  used  for  the  nozzle 
sizes  that  were  used  in  Figure  15  and  these  will  be  used  consistently 
throughout  the  remainder  of  the  text.  The  data  of  Figure  l6  axe  the 
same  data  that  are  plotted  in  Figure  15 .  A  comparison  of  the  two  figures 
indicates  the  advantage  of  using  the  dimensionless  coordinates  to  elimi¬ 
nate  the  pressure  effect.  The  data  of  Figures  l6,  17,  and  18  extend 
over  a  limited  range  of  frequencies  and  L*,  and,  therefore,  it  is 
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Figure  15.  Frequency  versus  l/L*  for  TF  propellant.  Lines  have  been  drawn 
through  data  taken  at  approximately  constant  pressure. 
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difficult  to  determine  the  exact  shape  of  the  curve.  However,  in  the 
following  chapter  the  general  nature  of  the  curve  will  he  postulated. 

It  should  he  noted  that  small  values  of  l/£  actually  represent 
large  values  of  L*  and  are  therefore  indicative  of  a  more  stable 
regime.  Because  of  this,  there  are  very  few  data  points  for  l/l  values 
of  about  30  or  less.  As  the  value  of  l/l  becomes  larger,  the  value  of 
L*  becomes  smaller  approaching  the  limiting  conditions  of  the  low-L* 
burner.  In  terms  of  experimental  results,  this  corresponds  to  a  region 
of  erratic  chuffing  where  it  was  difficult  if  not  impossible  to  obtain 
firings  containing  oscillations  that  were  consistent  enough  to  analyze 
in  terms  of  the  frequency  and  L*. 

Figure  19  is  a  plot  of  G  propellant  data  on  7  -  l/l  coordinates. 

It  is  readily  observed  from  this  figure  that  there  is  a  pressure  effect 
of  some  nature.  However,  the  higher  frequency  data  represent  three 
different  nozzles  and  a  pressure  range  of  between  approximately  20  to 
40  psia,  whereas  the  lower  curve  (excluding  the  data  of  GB  propellant) 
is  data  representing  one  nozzle  size  only,  corresponding  to  a  pressure 
of  20  psia.  A  precise  explanation  as  to  the  radical  change  in  frequency 
with  the  change  in  pressure  is  not  readily  available  at  this  time. 

However,  as  was  discussed  briefly  in  Section  C  of  this  chapter,  data 
obtained  using  GB  propellant  indicate  that  the  transmission  of  radiant 
energy  through  G  propellant  probably  is  an  important  factor.  Two  of 
the  data  points  obtained  with  GB  propellant  represent  firings  made  at 
45  to  50  psia,  and  even  at  these  comparatively  high  pressures  the  fre¬ 
quencies  for  the  run  were  in  the  lower  range,  thus  implying  that  radiation 
does  cause  the  large  increase  in  the  frequency. 
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In  order  to  shov  the  relative  location  of  the  different  data  from 
the  various  propellants,  the  data  of  TF,  F,  XF,  and  GB  propellants 
have  been  plotted  together  on  semi -logarithmic  coordinates  in  Figure  20 
(l/6  being  the  logarithmic  scale  in  order  to  better  include  the  data 
of  GB  propellant) .  The  data  of  the  three  catalyzed  propellants  seem 
to  fall  on  essentially  the  same  curve  with  the  GB  propellant  data  at 
the  extreme  of  the  curve.  The  data  of  the  least  stable  of  the  propel¬ 
lants,  XF,  are  found  at  the  left  of  the  plot  with  the  GB  propellant 
data  representing  the  most  stable  propellant  on  the  right  of  the  plot. 

It  is  impossible  to  compare  the  results  of  the  frequency-L* 
plots  with  results  of  other  investigators  because  data  of  this  type 
have  not  been  reported  before.  Although  other  investigators  have 
reported  a  qualitative  influence  of  L*  in  studies  of  NAI,  prior  to 
the  present  investigation  a  quantitative  dependence  of  non-acoustic 
frequencies  on  L*  has  not  been  reported.  The  followxi.g  chapter  contains 
a  comparison  of  the  results  of  the  theoretical  analysis  made  in  Chapter 
II  with  the  experimental  results  that  have  been  presented  in  this 


chapter . 
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CHAPTER  V 


A  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
FOR  TF  PROPELLANT 

This  chapter  is  devoted  to  a  discussion  of  the  results  obtained 
from  a  parametric  study  of  Equations  (2.19)  and  (2.27),  and  also  an 
interpretation  of  the  observations.  The  data  obtained  from  TF  pro¬ 
pellant,  the  most  thoroughly  studied  of  the  five  investigated,  have 
been  used  as  a  reference  by  which  the  results  of  the  computations 
are  Judged. 

It  will  be  recalled  from  Chapter  II  that  Equations  (2.19)  and 
(2.27)  are  analogous  equations  representing,  respectively,  the  cases 
for  which  kinetic  factors  and  thermodynamic  factors  regulate  the  sur¬ 
face  reactions.  Both  equations  are  of  the  form 

l/5  »  7f  ( 7 ,  q,  parameters)  (5.l) 

where  f  (7,  t),  parameters)  represents  two  different  functions  as 
defined  by  Equations  (2.19)  and  (2.27).  These  equations  formed  the 
basis  of  the  parametric  study  and  were  programed  to  be  rua  on  the  IBM 
7040  digital  computer  of  the  University  of  Utah  Computer  Center. 

The  fortran  statements  for  the  two  programs  have  been  recorded  as 
Appendix  H  in  the  present  work. 
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A.  A  Discussion  of  the  Parameters 

The  parameters  that  were  varied  in  this  study  were  the  activation 
energy  (E),  the  mean  surface  temperature  (Tg),  and  the  heat  of  gasification 
(q).  It  has  been  noted  previously  that  activation  energies  have  been 
reported  everywhere  from  l6  to  50  Kcal/mole,  with  Nachbar  and  Williams 
[22]  calculating  a  value  of  approximately  30  Kcal/mole  from  a  mathe¬ 
matical  model  for  a  hot  plate  experiment  taking  into  account  a  possible 
discrepancy  in  the  surface  temperature  of  the  plate. 

Pawling  and  Smith  [23,24]  have  studied  the  surface  temperature 
of  burning  strands  of  ammonium  perchlorate  as  well  as  composite 
propellants.  Their  studies  indicate  that  the  surface  temperature  of 
the  burning  propellant  should  be  between  775  and  900°K  for  pressures 
ranging  from  20  to  200  psia. 

The  heat  of  gasification  at  the  burning  surface  is  not  as  well 
defined  as  the  other  terms  and  is  therefore  much  more  difficult  to 
measure  or  evaluate.  However,  a  qualitative  range  of  bounding  values 
can  be  approximated.  The  heat  of  sublimation  of  ammonium  perchlorate 
is  5 6  Kcal/mole  or  approximately  480  cal/gm  and  can  be  considered  as 
an  upper  limit  for  the  heat  of  gasification.  The  reactions  occurring 
at  the  surface  Involve  both  the  ammonium  perchlorate  and  the  fuel 
binder  reacting  together  as  well  as  the  thermal  decomposition  of 
the  ammo nl urn  perchlorate.  Heterogeneous  reactions  between  the  oxidizer 
and  fuel  would  tend  to  be  exothermic,  thus  causing  the  overall  heat  of 
gasification  to  be  considerably  smaller  (more  exothermic)  than  the 


heat  of  sublimation. 
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These  three  parameters  are  implicitly  related  to  and  connected 
with  the  surface  of  the  burning  propellant.  Therefore,  it  might  be 
worthwhile  to  interject  considerations  of  the  real  character  of  the 
propellant  surface.  While  it  is  difficult  to  determine  the  exact 
nature  of  the  surface  during  burning,  it  is  possible  to  interrupt 
the  burning  and  observe  the  quenched  surface. 

Figure  21  is  a  photomicrograph  of  the  surface  of  a  sample  of 
propellant  that  was  extinguished  during  chuffing.  The  large  solid 
crystals  of  the  perchlorate  can  be  seen  protruding  well  above  the 
fuel  binder  and  rest  of  the  propellant,  demonstrating  the  lack  of 
homogeneity  and  uniformity  at  the  surface. 

In  the  propellants  used  in  this  study,  fifty  per  cent  by  weight 
of  the  ammonium  perchlorate  in  the  propellant  was  coarse  material, 
crystals  ranging  between  150  to  300  microns  in  size.  Even  for  the 
slowest  burning  propellant,  the  heated  zone  adjacent  to  the  burning 
surface  extends  only  about  200  microns  into  the  solid.  Crystals, 
such  as  seen  in  Figure  21,  are  probably  at  the  initial  temperature  at 
one  extreme  and  protruding  well  into  the  flame  zone  at  the  other.  To 
discuss  a  value  for  a  surface  temperature  in  the  light  of  such  knowledge 
might  appear,  and  indeed  is,  naive.  However,  the  tractnbility  of  the 
problem  depends  upon  the  assumption  that  to  talk  about  a  surface  tempera¬ 
ture  is  meaningful.  Some  comfort  can  be  taken  from  the  fact  that  the 
concept  of  an  ignition  temperature,  to  which  the  same  objections  can 
be  raised,  has  proved  to  be  highly  successful  in  describing  the  ignition 
process  [for  example  see  reference  32]. 


Scale  — **| 
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Figure  21.  The  nature  of  the  burning  surface.  A  photomicro¬ 
graph  of  the  surface  of  a  sample  of  F  propellant  which  extinguished 
after  burning  unstably. 
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The  difficulty  in  evaluating  the  heat  of  gasification  can  now 
be  tetter  understood.  The  heat  of  gasification  is  the  sum  of  heats  of 
reaction  and  sublimation  for  all  reactions  in  vhich  a  solid  phase 
participates,  including  pyrolysis  of  the  binder,  the  sublimation  or 
thermal  decomposition  of  the  perchlorate  crystals  themselves,  and  what¬ 
ever  heterogeneous  reactions  occur. 

In  the  following  sections  of  this  chapter,  variation  of  these 
parameters  to  produce  agreement  between  Equations  (2.19)  and  (2.27) 
and  the  data  will  be  discussed. 

B.  Theoretical  Results  for  the  Kinetic  Limited  Case 

It  has  been  noted  severed,  times  that  Equation  (2.19)  ( or  (5-l)) 
describes  l/6  as  a  function  of  7,  n,  and  the  parameters.  The  assumption 
that  the  surface  temperature  is  determined  by  kinetic  factors  introduces 
an  activation  energy  into  the  theory  (see  Equation  (2.2l)),  and  this 
activation  energy  will  therefore  be  included  in  the  discussion  of  the 
parameters .  Because  the  stability  limit  is  of  primary  concern,  tj  has 
been  taken  as  equal  to  zero  in  the  discussion  that  follows  except 
where  otherwise  noted.  The  phase  shift  between  the  pressure  and  the 
energy  flux  that  appears  as  tangent  5  in  the  equations  has  also  been 
taken  as  zero  except  where  noted. 

After  varying  the  value  of  the  heat  of  gasification  over  a 
range  of  values  between  -100  to  400  cal/gm,  it  was  determined  that 
values  of  q  near  zero  gave  the  best  fit  of  the  computer  solution  with 
the  experimental  results .  Figure  22  demonstrates  the  dependency  of 


the  theoretical  curve  on  the  value  of  q  when  other  parameters  are  held 
constant. 

One  should  keep  in  mind  that  the  stable  regime  for  combustion 
lies  to  the  left  of  the  paper,  corresponding  to  large  values  of  L*, 
and  with  the  growth  factor  taken  as  zero,  the  curves  in  Figure  22 
correspond  to  stability  boundaries.  The  point  where  the  actual  stability 
boundary  passes  through  the  data  points  can  be  determined  from  the 
experimental  results.  Two  of  the  triangular  points  in  Figure  22 
corresponding  to  the  approximate  coordinates  of  7  equal  11  and  l/l  equal 
30  were  actually  located  in  the  stable  regime  and  have  negative  growth 
factors .  The  remainder  of  the  points  on  the  plot  correspond  to  data 
obtained  in  the  unstable  regiu.e  with  a  few  neutral  points .  Considering 
this  and  utilizing  the  results  of  Figure  13  (the  correlation  of  growth 
factor  data),  the  stability  limit  can  be  approximately  located  at  the 
coordinates  mentioned  above  (i.e.,  7  =  11,  l/£  =  30 ). 

With  this  information  and  referring  again  to  Figure  22  it 
appears  that  the  value  of  q  which  agrees  best  with  the  experimental 
results  is  be  .ween  0  and  20  cal/gm.  (depending  of  course  on  the  values 
of  E  and  Tg).  The  value  for  the  activation  energy  was  varied  between 
20,000  and  40,000  kcal/mole  with  the  results  represented  in  Figure  23. 

The  best  fit  with  the  experimental  results  yields  a  number  between 
30,000  and  40, 000  Kcal/mole.  The  former  may  be  a  more  reasonable  value 
to  select  considering  that  it  coincides  with  that  reported  by  Nachbar 
and  Williams  [22  ] . 

Figure  24  demonstrates  the  effect  of  varying  the  mean  surface 
temperature  between  750  and  950  K.  From  the  results  of  this  figure 
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it  can  be  concluded,  that  the  solution  of  Equation  (2.19)  is  only  slightly 
sensitive  to  variations  in  the  surface  temperature.  An  average  value 
of  85O  K  was  therefore  used  in  the  remainder  of  the  calculations  without 
taking  into  account  the  small  variation  in  temperature  that  occurs 
with  change  in  pressure  level. 

The  computer  curves  of  Figures  22,  23>  and.  2  k  correspond  to 
solutions  with  the  growth  factor  taken  as  zero.  In  the  previous  dis¬ 
cussion  it  was  pointed  out  that  the  growth  term  is  zero  only  at  one 
point  of  the  dimensionless  L*-frequency  domain,  the  stability  limit. 

In  order  to  determine  the  effect  of  q  on  the  theoretical  results,  an 
empirical  relationship  between  q,  7  and  I  was  obtained  by  means  of 
Figure  13  in  the  preceding  chapter.  If  the  straight  line  relationship 
that  was  discussed  is  considered,  then  Equation  (5-2)  follows  as  the 
mathematical  expression  corresponding  to  Figure  13 


This  equation  has  been  plotted  in  Figure  25  along  with  the  computer 
solution  to  Equation  (2.19)  using  what  appears  to  be  the  best  values 
of  the  parameters.  The  intersection  of  these  two  curves  should  define 
the  stability  limit.  If  the  value  of  q  is  taken  as  zero,  the  inter¬ 
section  of  the  computer  solution  with  Equation  (5*2)  occurs  at  approxi 
mately  7  =  10  and  l/i  =  25  instead  of  the  location  observed  in  Figure 
25.  From  the  above  discussion  of  the  location  of  the  stability  limit, 
it  can  be  seen  that  the  value  of  q  should  be  between  0  and  20  cal/gm. 
(assuming  that  the  other  parameters  are  correct).  The  exact  location 


Figure  25.  Comparison  oi‘  theoretical  results  (tj=0)  with  the  empirical 
results  obtained  from  growth  factor  data  {the  parabolic  curve). 
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of  this  intersection  is  determined  by  Equation  (5*2)>  an  empirical 
relationship,  because  a  theoretics!,  basis  for  locating  it  has  not  yet 
been  achieved. 

It  will  be  recalled  that  Equation  (2.19)  includes  a  phase  shift 
between  the  pressure  and  the  energy  flux  which  was  taken  as  zero  in 
Figures  22-25.  The  results  of  taking  the  tangent  of  the  phase  angle 
as  0.3,  and  the  assumed  best  val-ues  of  the  activation  energy  and  the 
surface  temperature  from  the  parametric  study,  have  been  plotted  as 
Figure  26.  The  shape  of  the  curves  for  various  values  of  q  are  similar 
to  those  where  the  phase  shift  is  taken  as  zero,  but  larger  values  of 
q  are  now  needed  to  fit  the  experimental  results  better.  A  small 
negative  phase  shift  was  also  considered,  but  there  were  no  solutions 
to  the  equation  for  these  conditions.  These  results  indicate  that  in 
order  to  obtain  meaningful  results  from  Equation  (2.19),  either  the 
energy  flux  must  lead  the  pressure  by  a  small  positive  phase  angle 
or  the  two  are  in  phase.  This  is  in  agreement  with  the  discussion  in 
Chapter  II,  where  the  relationship  between  the  flux  and  the  pressure 
was  considered  in  some  detail.  It  was  concluded  there  that  a  phase 
shift  should  be  small  or  zero. 

C.  Theoretical  Results  for  the  Equilibrium  Vaporization  Case 

The  assumption  that  the  surface  temperature  is  determined  by 
equilibrium  vaporization  and  subject  to  evaluation  by  the  Clausius- 
Clapeyron  equation,  introduces  a  heat  of  sublimation  into  the  analysis 
(see  Equations  (2.2U)  and  (2.27)).  The  value  of  the  heat  of  sublimation 
for  ammonium  perchlorate  has  been  reported  in  the  literature  as  being 


ooo'oe  = 


111 

(V  r-  lO  O  00 

1  l  1 

_l  — . 

o  oo  io  in  cm 

N  Z 

CO  CM  CM  (\|  CO 

N  3 

o 

o  o  O  O  O 

z 

_J 

§ 

4  O  O  ■  • 

> 

co 

o 

in  3 

qo 

II  c 

m  .A) 


Figure  26.  Effect  of  including  a  phase  shift  on  theoretical  results  (tj=0) 
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56  Kcal/mole  [22,23,2k],  and  will  not  be  subjected  to  the  parametric 
study  but  will  be  considered  as  a  known  constant.  Powling  and  Smith, 
who  are  among  those  that  propose  the  thermodynamic  mechanism,  report 
surface  temperatures  as  a  function  of  pressure  [24].  They  report  values 
between  775  and  900  K  for  pressures  up  to  200  psia,  and  the  appropriate 
values  were  used  in  this  particular  study,  thus  excluding  another 
parameter,  and  leaving  only  the  heat  of  gasification  to  be  varied. 

Figure  27  contains  the  results  of  the  analysis  carried  out 
varying  q  in  Equation  (2.27).  The  curves  for  q  equal  to  80  and  100 
cal/gm.  appear  to  pass  nearest  the  location  of  the  stability  point. 
Utilization  of  Equation  (5.2)  with  the  results  of  Figure  27  indicates 
that  the  assumption  of  equilibrium  vaporization  appears  to  be  valid, 
at  least  for  TF  propellant. 

The  results  of  this  chapter  have  been  for  TF  propellant  only. 

The  extension  and  applicability  of  these  results  to  the  other  pro¬ 
pellants  used  in  this  study  and  other  classes  of  propellants  in  general 
will  be  considered  in  the  next  chapter. 


CHAPTER  VI 


THE  OVERALL  SCOPE  OF  THE  THEORETICAL  RESULTS 

The  parameters  of  the  theoretical  equations  have  been  adjusted 
to  fit  the  data  for  TF  propellant.  We  inquire  now  as  to  the  appli¬ 
cability  of  these  parameter  values  to  the  other  propellants  of  this 
study  or  to  other  classes  of  propellants .  An  examination  of  Equations 
(5.l),  (2.19),  and  (2.18)  reveals  that  the  parameters  involved  are 
the  activation  energy,  the  mean  surface  temperature,  the  initial  tempera¬ 
ture  of  the  propellant,  the  heat  of  gasification,  and  the  heat  capacity 
of  the  solid.  The  initial  temperature  iB  arbitrary  and  can  be  con¬ 
sidered  as  a  constant  for  present  purposes.  The  heat  capacity  varies 
slightly  among  propellants  and  actually  enters  into  the  calculations 
in  a  rather  minor  way,  so  that  variations  in  this  quantity  can  be 
ignored.  This  leaves  the  three  parameters  that  have  already  been  dis¬ 
cussed  in  some  detail.  It  can  therefore  be  concluded  that  the  curves 
in  Figures  22  through  27  apply  to  any  propellant  whose  parameters  have 
the  values  specified  for  the  individual  cases.  The  questions  now 
arise  as  to  whether  these  parameter  values  vary  from  propellant  to 
propellant,  and  if  so,  to  what  extent. 

A.  The  PBAA  PropeJlants  Used  in  this  Study 

For  convenience,  the  definition  of  the  dimensionless  L*  and 
frequency  are  repeated  here  as  Equation  (6.l) 
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,  C*r  ~L» 

5  “  k*RTf 


(6.1) 


From  these  definitions  it  can  be  seen  that  other  properties  of  the 
propellant,  particularly  the  burning  rate,  enter  into  the  calculations 
of  the  dimensionless  variables  and  therefore  contribute  to  the  location 
of  the  points  along  the  curve. 

If  it  is  assumed  that  propellants  F,  XF,  and  GB  have  the  same 
activation  energy,  surface  temperature,  and  heat  of  gasification  as 
propellant  TF  the  remarks  of  the  previous  chapter  apply  to  these  pro¬ 
pellants  also.  From  Figure  28  it  can  be  seen  that  the  results  for 
various  propellants  follow  the  same  trend,  the  difference  in  propellants 
being  mainly  in  the  general  location  of  data  points  along  a  common  cor¬ 
relating  line.  The  aluminized  propellants  fall  further  to  the  left  on 
the  dimensionless  plot,  while  the  non-alumlnlzed  propellants  are  shifted 
to  the  right.  It  will  be  recalled  that  the  termination  of  these  points 
determines  the  stability  limit,  which  is  of  prime  interest  and  will 
now  be  considered  further. 

In  Figure  28  the  upper  curve  (marked  q=20)  is  the  computer  solution 
of  Equation  (2.19)  with  an  activation  energy  of  30,000  Kcal/mole,  a 
surface  temperature  of  850  K,  a  heat  of  gasification  of  20  cal/gm., 
and  the  growth  factor  taken  as  zero.  We  assume  for  the  moment  that 
the  same  parameter  values  apply  to  F,  XF,  and  GB  propellants  as  well 
as  TF  propellant.  The  second  curve  is  the  empirically  determined  relation¬ 
ship  for  TF  propellant  that  was  obtained  from  the  growth  factor  data  for 
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Figure  28.  A  comparison  of  the_data  from  the  four  PBAA  propellants  with 
the  computer  solution  for  E  =  30,000,  T  =  850,  q  =  20,  and  q  =  0.  The 
empirical  relationship  of  Figure  13  has  also  been  included. 
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that  propellant.  Equation  (5.2).  The  Intersection  of  these  two  curves 
should  be  the  stability  limit  for  TF  propellant.  The  empirical  relation¬ 
ship  appears  to  follow  the  trend  of  all  the  data  quite  well,  so  that 
if  it  were  assumed  that  the  same  relationship  held  for  all  four  pro¬ 
pellants,  then  there  would  be  a  unique  stability  limit  for  the  four 
propellants.  The  experimentally  determined  fact  that  the  stability 
limit  varies  from  propellant  to  propellant  indicates  that  one  of  two 
possibilities  exists:  either  (l)  the  empirical  relationship  does  not 
hold  for  all  of  the  propellants,  or  (2)  the  values  of  the  parameters 
change  from  one  propellant  to  the  next.  Actually  a  third  possibility 
exists  in  that  both  of  the  preceding  two  could  be  true. 

Recalling  the  nature  of  the  heat  of  gasification  term,  it  could 
well  be  suspected  of  being  different  for  various  propellants .  The 
activation  energy  is  probably  different  also,  particularly  between 
catalyzed  and  uncatalyzed  propellants . 

The  data  at  hand  do  not  allow  obtaining  growth  factor  values 
for  propellants  other  than  TF.  The  fact  that  Equation  (5.2)  fits  the 
7  vs.  l/i  data  for  the  other  propellants  quite  well  could  inqsly  that 
an  expression  of  this  type  might  be  general  for  a  class  of  propellants. 

Further  investigation  is  necessary  to  determine  the  effect  of 
propellant  composition  upon  the  above-mentioned  parameters,  and  the 
nature  of  the  relationship  between  7,  t,  and  the  growth  factor. 

B.  Extension  of  the  Theory  to  Other  Propellants 

In  the  preceding  chapter  it  was  shown  that  variations  in  the 
surface  temperature  and  the  activation  energy  have  a  very  slight  effect 
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on  the  theoretical  results.  Also,  it  is  probably  true  that  the  values 
of  these  parameters  depend  to  a  large  extent  on  the  nature  of  the  oxidizer 
and  are  fairly  independent  of  the  binder  in  the  propellant  (for  example, 
see  [4]  and  [22]).  This  would  tend  to  indicate  that  the  values  used 
for  these  parameters  would  be  approximately  the  same  for  ammonium  per¬ 
chlorate  propellants  using  a  binder  other  than  PBAA.  Therefore,  the 
location  of  the  theoretical  curve  appears  to  be  dependent  on  the  value 
of  the  heat  of  gasification  more  than  any  other  parameter.  This  con¬ 
clusion  does  not  seem  unreasonable  when  one  considers  the  nature  of 
the  term.  It  must  be  dependent  upon  the  binder  of  the  propellant 
and  the  manner  in  which  it  reacts  with  the  oxidizer,  and  is  probably 
influenced  by  metal  additives  in  the  propellant  also. 

The  present  study  has  been  involved  with  composite  propellants, 
but  it  is  reasonable  to  assume  that  the  analysis  could  be  applied  to 
double-base  propellants  without  difficulty.  The  values  of  the  surface 
temperatures,  activation  energies,  and  heats  of  gasification  for  these 
propellants  are  considerably  different  from  those  for  composite  pro¬ 
pellants,  but  so  are  the  stability  limits  (L*  as  high  as  5000  cm.  [20]). 
The  lose  of  the  correct  values  for  these  parameters  with  Equation  (2.19) 
(evaluating  eQ/p  properly)  should  give  reasonable  results. 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

The  objectives  of  the  present  study  were  twofold:  to  investigate 
non-acoustic  instability  experimentally,  determining  stability  limits 
for  the  reference  propellants,  and  to  attempt  a  mathematical  explanation 
of  the  experimental  observations.  The  present  chapter  will  include  a 
discussion  of  these  two  objectives  plus  a  comparison  of  the  theoretical 
implications  of  the  study  with  the  experimental  observations  of  other 
investigators . 

A.  Experimental  Observations 

The  experimentally  determined  stability  limits  have  been  presented 
as  Figure  10.  The  results  of  this  figure  indicate  that  the  more  highly 
aluminized  propellants  are  more  unstable  and  the  uncatalyzed  G  and  GB 
propellants  appear  to  be  the  most  stable  of  the  propellants  in  terms 
of  the  pressure-L*  domain. 

It  was  observed  that  non-acoustic  frequencies  are  dependent  on 
the  value  of  L*  and  can  be  correlated  by  plotting  frequency  versus  the 
reciprocal  of  L*  (see  Figure  15).  Using  the  dimensionless  coordinates 
of  7  and  l/l,  the  pressure  effect  of  the  f requency-L*  dependency  was 
eliminated.  Figure  20  shows  the  data  obtained  vising  TF,  F,  XF,  and  GB 
propellants  plotted  on  the  dimensionless  coordinates.  The  data  from  the 
four  propellants  fall  within  approximately  the  same  bounds,  being  grouped 
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according  to  the  respective  propellants.  The  general  trend  is  for 
frequency  to  increase  as  L*  decreases. 

One  of  the  propellants  tested,  G,  behaved  quite  differently  from 
the  others  and  the  data  obtained  for  this  propellant  did  not  correlate 
with  those  obtained  for  the  others.  The  results  for  G  propellant  are 
graphically  displayed  in  Figure  19-  It  is  thought  that  this  discrepancy 
is  probably  due  to  the  effects  of  radiation  (G  propellant  is  slightly 
translucent);  therefore,  its  composition  was  modified  slightly  to 
include  carbon  black  in  an  attempt  to  eliminate  these  effects.  The 
data  of  GB  propellant,  that  containing  the  carbon  black,  do  correlate 
with  the  data  of  the  other  propellants,  and  therefore  it  is  felt  that 
the  discrepancies  observed  in  the  data  of  G  propellant  were  in  fact 
due  to  radiation. 

Assuming  that  the  growth  of  the  pressure  oscillations  was 
exponential.,  growth  factors  were  calculated  from  the  pressure  time 
traces .  A  correlation  between  these  and  the  dimensionless  frequency 
and  L*  indicated  a  relationship  almost  independent  of  the  growth  factor 
(see  Equation  (5*2)  and  Figure  13).  This  correlation  was  used  to  lo¬ 
cate  the  stability  limit  for  TF  propellant  with  respect  to  the  dimension¬ 
less  coordinates. 

B.  Theoretical  Results 

The  theoretical  analysis  that  was  developed  in  Chapter  II  is  based 
on  an  energy  balance  for  the  solid  propellant  considering  energy  accumu¬ 
lation  in  the  solid  and  a  mass  balance  on  the  rocket  chamber  considering 
mass  accumulation  in  the  gas  phase.  The  analysis  is  based  on  the 
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assumption  that  oscillations  grow  from  small  perturbations.  Two  other 
assumptions  which  are  very  important  in  the  analysis  are:  (l)  the 
burning  surface  is  homogeneous  and  uniform,  and  (2)  the  phase  lag 
between  the  pressure  and  the  energy  flux  is  small  or  zero. 

The  analysis  thus  developed,  reduces  to  e~i  expression  of  the 

form 

l/fi  =  7  f  (7,  »b  parameters)  (5-l) 

The  principle  parameters  involved  are  the  activation  energy,  the  mean 
surface  temperature,  and  the  heat  of  gasification.  The  parametric 
study  discussed  in  Chapter  V  indicates  that  an  activation  energy  of 
30,000  Kcal/mole,  a  surface  temperature  of  850  K,and  a  heat  of  gasifi¬ 
cation  of  between  O  and  20  cal/gm.  result  in  the  best  fit  of  the 
experimental  results,  although  other  combinations  of  the  parameter 
values  fit  quite  well  also.  The  values  of  the  first  two  of  these 
parameters  are  conq) arable  to  those  reported  in  the  literature.  The 
value  for  the  heat  of  gasification  appears  to  be  the  right  order  of 
magnitude  even  though  this  quantity  has  not  been  measured  experimentally. 

The  fact  that  the  theoretical  analysis  results  in  an  expression 
that  agrees  quantitatively  with  the  experimental  results  tends  to 
vindicate  the  analysis  and  its  attendent  assumptions. 

The  analysis  was  developed  considering  composite  propellants 
but  should  be  applicable  to  double-base  propellants  as  well. 

C .  A  Comparison  with  the  Results  of  Other  Investigators 

The  theoretical  results  of  this  study  with  the  best-fit  para¬ 
meters  are  concisely  presented  in  Figures  25  and  28,  and  can  be  used 
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to  describe  the  results  reported  by  others  in  a  qualitative  manner . 

Considering  first  the  theory  developed  by  Akiba  and  Tanno  and 
extended  by  Sehgal  and  Strand,  it  is  observed  that  with  the  aid  of  the 
empirical  relationship  relating  the  growth  factor  to  the  frequency 
and  L»,  Figure  13,  %  stability  limit  can  be  determined.  The  value 
of  L*  at  that  limit  is 


L»  =  - £ 

cr  C*  ?2 


(7.1) 


This  is  essentially  the  same  equation  that  is  developed  by  the  above 

authors  with  i  replaced  by  a  theoretically  determined  x  .  Thus 
cr  cr 

the  present  analysis  appears  to  be  compatible  with  the  results  of  the 
Akiba-Tanno  theory. 

Yount  and  Angelus  [40]  have  observed  little  or  no  dependency  of 
frequency  upon  L*.  It  would  appear  to  the  author  that  they  are  probably 
operating  at  small  i  values  where  it  appears  (see  Figure  25)  that  the 
frequency  exhibits  very  litttle  dependency  on  L*. 

Price  et  al.  report  data  from  their  lov-L*  burner  in  terms  of 
the  preferred  frequency  concept  that  hew  been  discussed  to  some  extent 
already  (see  Figure  14).  This  also  is  in  qualitative  agreement  with 
the  results  of  the  present  investigation.  Figure  15  is  a  plot  of 
frequency  versus  the  reciprocal  of  L*  and  the  lines  through  the  data 
represent  constant  pressure  lines  for  an  idealized  system  (as  predicted 
by  the  theory).  If  one  follows  along  a  constant  pressure  curve  beginning 
at  the  stability  limit  (which  corresponds  to  a  particular  frequency), 
the  frequency  increases  as  L*  decreases  until  what  hew  an  appearance  of 


-  95  - 


f 


being  an  asymptotic  limit  is  reached  for  a  very  small  L* .  One  obtains 
a  band  of  frequencies  having  values  between  those  corresponding  to  the 
stability  limit  and  the  "asymptotic"  limit.  An  increase  in  the  pres¬ 
sure  would  result  in  another  band  of  frequencies  at  a  higher  level, 
the  result  being  the  preferred  frequency  concept. 

In  conclusion  it  can  be  noted  that  although  the  observation  and 
prediction  of  a  frequency-L*  dependency  reported  in  the  present  investi¬ 
gation  may  appear  somewhat  radical  (in  that  it  has  not  been  reported 
before),  it  is  not  contradictory  with  the  results  of  previously  reported 
work.  In  fact  it  is  complementary  to  the  results  of  other  investigators 
and  tends  to  correlate  somewhat  the  apparently  unrelated  results  that 
have  been  reported  up  to  the  present  time. 
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APPENDIX  A 


NOMENCLATURE 

Dimensions:  L  =  length,  M  =  mass,  0  =  time,  and  T  =  temperature 


Roman  Letters 

At 

A 

C* 


E 

f 

F 


H  + 
s  + 


H 

i 

Im(  ) 


area  of  burning  surface 

nozzle  throat  area 

Arrhenius  frequency  factor 

characteristic  velocity 

nozzle  coefficient 

specific  heat  at  constant  pressure 

activation  energy 

incident  energy  flux 

m 


M 

Im( 


c  T  9. 
p  a  0 

HP 

?cp~Vn 

2HSP 

temperature  gradient  at  the 
surface 

temperature  gradient  on  the  gas 
phase  side  of  the  surface 

«  +  cp(fs  "  T0) 


(-Dl/2 


L 

L/e 

L/0 

1/0 

L 2/T102 

L2/02 

M/03 

dimensionless 

dimensionless 

dimensionless 


T/L 


t/l 

l2/92 


dimensionless 


imaginary  part  of  a  complex 
expression 


101  - 


102 


Roman  Letters  (cont'dj 


ratio  of  the  burning  surface  area 
to  the  nozzle  throat  area 

thermal  conductivity 


dimensionless 


ML/T9- 


time  constant 

ratio  of  motor  free  volume  to 
nozzle  throat  area 


molecular  weight 


dimensionless 


burning  rate  exponent 


dimensionless 


pressure 


ambient  pressure  into  which 
the  rocket  exhausts 

heat  of  phase  change 

heat  of  gasification 


L2/e2 

L2/e2 


burning  rate 
gas  constant 


l2/92t 


Re(  ) 


Real  part  of  a  complex 
expression 

,  ,  \2  2 
(a  +  l)  +  t 


dimensionless 


temperature 

adiabatic  flame  temperature 


temperature  at  x  =  °° 


volume 


ni  +  t)> 


dimensionless 


distance  into  the  propellant 

2c  (T  -  T  ) 
p  s  o 


dimensionless 
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Roman  Letters 


Y 


1 


(cont'd) 

7c  RT2 
-P— g 
2HSAH 

s 


c  RT2 
P  5 
HAH 

s 


Greek  Letters 


O 

a 

g 

0 


7 

r 

AH 

€ 

s 


e 


1 

e 

9< 

x 


thermal  diffusivity 
growth  factor 

burning  rate  perturbation 

hctm 

-2 

r 

ratio  of  specific  heats 

heat  of  sublimation 

volume  fraction 

phase  shift  between  flux 
and  pressure 

growth  factor  (ag/o>) 

temperature  perturbation 

surface  temperature  perturbation 
r 

2a 


X) 

t 

K 

P 

a 


ratio  of  thermal  conductivities 

-2 

C*r 

L*,  nondimensionalized  by  •j— ^ 
pressure  perturbation 
density 

dumny  variable  defined  by 
Equation  (2.11c) 

dummy  variable  defined  by 
Equation  (2. lid) 


dimensionless 


dimensionless 


L2/0 

1/0 

dimensionless 

dimensionless 

dimensionless 

L2/02 

dimensionless 

radians 

dimensionless 

dimensionless 

dimensionless 

l/L 

dimensionless 

dimensionless 

dimensionless 

m/l3 

dimensionless 

dimensionless 


Greek  Letters  (cont'd) 


4>  energy  flux  perturbation  dimensionless 

a)  angular  frequency  l/6 


Subscripts 

cr 

eh,  g 

P 

s 

t 


refers  to  the  stability  limit 

refers  to  a  condition  or  property  in  the 
gas  chamber 

refers  to  a  property  of  the  propellant 
refers  to  the  gas-solid  interface 
refers  to  conditions  at  the  nozzle  throat 


Superscripts 

Bar  over  a  symbol  (e.g.  r)  refers  to  the  time-averaged  value 
Note:  Additional  nomenclature  specific  to  Appendix  G  has  been  listed 


in  that  Appendix. 


APPENDIX  B 


PROPELLANT  PROPERTIES  AND  COMPOSITIONS 

Propellant  Compositions 

PBAA  (polybutadiene- acrylic  acid  polymer)  and  ammonium  perchlorate 
are  the  basic  constituents  of  the  propellants  considered  in  the  present 
investigation.  The  propellants  each  contain  82  per  cent  (by  weight)  dispersed 
solids  and  18  per  cent  polymeric  binder;  however,  several  variations 
were  obtained  through  the  use  of  various  additives.  The  enact  compo¬ 
sition,  as  well  as  cor  •ooporli  ig  designation  of  each  propellant,  is  recorded 
in  Table  1.  The  procedure  and  conditions  for  compounding  the  propel¬ 
lant  are  as  follows :  the  propellant  was  mixed  at  6 5°C  at  an  absolute 
pressure  of  7  to  8  mm.  Hg  for  60  minutes,  cast  at  a  pressure  of 
7  to  8  mm.  Hg  and  cured  at  80°  for  7  days. 

The  composition  of  the  polyurethane  propellant  that  was  used  was 
obtained  from  the  Naval  Ordnance  Test  Station  [31]  and  is  recorded  in 
Table  2.  The  composition  was  modified  from  that  obtained  from  NOTS  to 
the  extent  that  the  ratio  of  coarse  to  fine  ammonium  perchlorate  was 
changed  from  50:50  to  70:30  in  order  to  reduce  the  viscosity  of  the 
uncured  propellant  during  the  mixing  and  casting  process.  The  mixing 
procedure  that  was  used  is  as  follows:  the  estsne,  1, 4-BD,  TEA,  and 
coarse  NH^CIO^  were  blended  together  in  a  pre-mix  for  15  minutes;  the 
fine  NH^CIO^  was  then  added  and  mixed  for  10  minutes;  finally  the  TMP 
weis  elded  and  mixing  continued  for  45  minutes.  The  mixing  temperature 
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TABLE  1 


COMPOSITION  OF  PBAA  PROPELLANTS 


(a) 


Solid  Ingredients  (weight  per  cent) 


NH^ClOi^ 


Pronel- 

lant 

Code 

Coarse^ 

15  micron^0 ^ 

Copper 
Chromite.  . 
Catalyst^' 

Carbon 

Black(e' 

(f) 

Aluminum 

G 

4l 

4l 

C 

0 

0 

GB 

40 

40 

0 

2 

0 

F 

40 

40 

2 

0 

0 

TF 

37-5 

37-5 

2 

0 

5 

XF 

35 

35 

2 

0 

10 

(a)  All  of  the  propellants  in  this  table  had  a  PBAA  (Thiokol,  poly- 
butadiene-acryllc  acid  copolymer)  binder  composed  of  85  per  cent 
(by  weight)  PBAA  and  15  per  cent  epoxy  resin  (Shell  Chemical  Co., 
Epon  828 ) . 

(b)  American  Potash  and  Chemical  Corp.,  designated  as  -48  +100. 

(c)  American  Potash  and  Chemical  Corp.,  designated  as  50^  less  than 
15  micron. 

(d)  Harshaw  Chemical  Co.,  Cu-0202  P. 

(e)  Cabot  Corp.,  Sterling  VR  carbon  black 

(f)  Reynolds  Aluminum,  No.  1-5H  atomized  powder  (50  per  cent  less 
than  25  micron). 
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TABLE  2 

COMPOSITION  OF  POLYURETHANE  PROPELLANT 


Weight  per  cent 

Ingredient 

56 

(a) 

-1+8  +100  ammonium  perchlorate 

2k 

10  micron  ammonium  perchlorate^^ 

19.32 

( c ) 

Estane' 

0.46 

\(c) 

IMP  ( Trimethyl ol  propane) 

0.14 

1, 1+-HD  (Butanediol)^ 

0.08 

TEA  (Triethanolamine)^0^ 

(a)  See  table  1,  footnote  (b). 

(b)  American  Potash  and  Chemical  Co.,  designated  as  10  micron 
(50  weight  per  cent  less  than  10  micron) 

(c)  All  of  the  Ingredients  for  the  binder  were  received  from  the 
Naval  Ordnance  Test  Station  at  China  Lake,  California. 
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was  75°C  and  both  the  mixing  and  casting  pressure  were  8  mm.  Hg. 

The  propellant  was  cured  for  two  days  at  80°C. 

Many  thermal,  chemical,  and  physical  properties  of  the  propel¬ 
lants  were  needed  in  the  various  calculations  throughout  the  course 
of  the  study.  The  methods  of  measuring,  calculating,  or  approxi¬ 
mating  these  properties  will  be  discussed  in  the  following  paragraphs. 

Density 

The  density  of  the  propellants  was  determined  by  measuring  a 
cylindrical  piece  of  propellant  approximately  1.5  inches  long  by 
1.5  inches  in  diameter  with  a  micrometer  and  then  weighing  the  sample 
on  an  analytical  balance.  The  densities  of  the  cured  binders  were 
measured  by  water  displacement,  and  the  density  of  ammonium  perchlorate 
was  obtained  from  the  International  Critical  Tables. 

Heat  Capacity 

The  heat  capacity  of  the  cured  binders  and  of  F  and  G  propellant 
was  determined  with  a  Dewar-flask-calorimeter  as  per  Ryan,  et  al.  [32]. 
Those  of  TF  and  XF  propellants  were  calculated  from  the  known  heat 
capacities  and  compositions  of  their  constituents. 

Thermal  Diffusivity 

The  thermal  diffusivity  of  F  and  G  propellant,  the  cured  binders, 
and  ammonium  perchlorate  was  the  same  as  that  determined  by  Ryan  et  al. 
[32],  and  were  determined  by  a  transient  technique  dependent  upon 
the  temperature-time  relationship  at  the  center  of  a  sample.  The 
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thermal  diffusivity  of  XF  and  TF  propellant  was  calculated  from  the 
equation 

n  k 

Pc 

Thermal  Conductivity 

The  thermal  conductivity  of  TF  and  XF  propellant  was  calculated 
from  the  Maxwell  equation  for  the  electrical  resistivity  of  heterogeneous 
materials  as  modified  for  thermal  conductivities  and  reported  by  Gorring 
and  Churchill  (l4).  This  equation  is 


k_  =  2  +  u  -  2e  (l-u) 
k.,  2  +  u  +  e (l-u) 


(B-l) 


where  v>  is  the  ratio  of  the  thermal  conductivities  of  the  discontinuous 
phase  (the  aluminum)  to  that  of  the  continuous  phase  (F  propellant), 
e  is  the  volume  fraction  of  the  discontinuous  phase,  and  is  the 
thermal  conductivity  of  the  continuous  phase.  The  thermal  conductivity 
of  an  aluminized  propellant  whose  thermal  diffusivity  had  been  measured 
by  the  method  outlined  above  was  checked  by  use  of  Equation  (B-l) 
and  found  to  be  within  4  per  cent  of  the  measured  value,  indicating 
that  the  equation  is  probably  an  excellent  approximation. 

The  measurements  of  the  density,  the  heat  capacity,  and  the  thermal 
diffusivity  are  probably  all  within  2  per  cent.  The  error  in  the 
approximation  for  the  thermal  conductivity  from  Equation  (B-l)  should 
be  less  than  10  per  cent. 
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Properties  of  the  Combustion  Products 

The  adiabatic  flame  temperature,  the  characteristic  velocity, 
and  the  ratio  of  the  specific  heats  for  the  combustion  products  of 
F,  G,  and  TF  propellant  were  calculated  considering  thermodynamic 
equilibrium  of  these  products.  The  results  of  these  calculations  were 
obtained  for  the  University  of  Utah  by  Dr.  Ralph  Coates  [9]  of  the 
Lockheed  Propulsion  Co.  The  average  molecular  weight  of  the  combustion 
products  was  then  calculated  from  the  definition  of  the  characteristic 
velocity 


c*2  = 


RT. 


r+i 

“■ter' r_1 


(B-2) 


The  values  of  the  above  properties  have  been  tabulated  and  appear 
in  Table  3.  It  was  assumed  that  the  properties  of  GB  propellant  were 
essentially  the  same  as  for  G  propellant,  and  those  of  XF  propellant 
were  estimated  from  those  of  TF  propellant. 


Linear  Burning  Rate 

The  steady-state  burning  rate  of  the  propellant  enters  into  the 
equations  of  Chapter  II  in  several  places,  and  is  therefore  a  very 
important  property.  Data  on  the  steady  burning  rate  for  the  various 
propellants  was  obtained  from  a  constant  pressure  strand  bomb.  Pro¬ 
pellant  strands  l/4  x  l/U  inch  and  approximately  3  inches  long  were 
burned  in  an  atmosphere  of  nitrogen  to  provide  these  data.  The  results 
are  plotted  in  Figure  29  showing  the  relationship  between  the  steady 
burning  rate  and  pressure  for  G,  F,  TF,  and  XF  propellants.  Several 
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PROPERTIES  OF  THE  PROPE: 
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gB  _ 1^60 _ .311  2.65  1*810  2630  1.21* 

(a)  Calculated  from  the  specific  heats  of  the  ingredients. 

(b)  Calculated  from  Equations  (B-2)  and  (B-l). 

(c)  Approximated 

(d)  The  properties  of  GB  propellant  were  assumed  to  be  the  same  as  those  for  G. 


Rate  (in/ sec) 


-  112 


Figure  29.  Burning  rate  versus  pressure  for  F,  TF,  XF,  and  G 
propellants .  Data  of  GB  propellant  fall  on  the  same  line  as  for  G 
propellant. 
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data  points  were  obtained  for  GB  propellant  and  fall  essentially  on 
the  line  determined  for  G  propellant.  It  was  therefore  assumed  that 
their  burning  rates  are  the  same.  Raw  data  have  been  included  in 
Table  19,  Appendix  I. 


APPENDIX  C 


EVALUATION  OF  THE  INTEGRAL  /  9/P  dx 


In  perturbing  the  energy  balance,  the  integral  of  9/p  as  a 
function  of  x  evolved  in  the  form 

c  T  9r 


.  C  CD  T  /Oo  C  i  W- 


(2.6) 


The  perturbed,  unsteady  heat  conduction  equation  provided  an  expression 
for  9/p  as 


Q(*) 

P 


\  +  ^s  -  T0} 

p  r(i-m)  t 


e-X(l+a+it)x  TQ)  e"2Xx 

7  (i+n)  t. 


(2.1) 


’ 

®0  k^s  -  T0}1 

-X(l+o+ix)x  k(T0  -  Tq)  _2Xx  ^ 

P  7(i+n)  T 

_  o  _ 

7( i+njTg 

(c-l) 

r  dx 


For  the  sake  of  convenience,  let  W  a 


U(T  -  T  ) 

B  °  .  Equation  (C-l) 


7(i+n)  T 


reduces  to 


pv-<w; 


-X(c+l+it )x 


dx 


•  W  / 
1  n 


-2Xx  . 
e  dx 


(C-la) 


Carrying  out  the  integration  leaves 


-  ilk  - 
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®0  -x(cj+1+1t)x  1  °°  . 

t  0  .  r,\  e  W  -2Xx 

(  P+W)  -Ma+l+ix)  o  +2X  e  0 


(C-2) 


Ml+o+iTj 


3>.» 


(l+o-if)  -  ^ 


(C-3) 


where  S  E  (l+o)2  +  t2 


Substitution  of  (C-3)  into  (2.6)  gives 


p  “  1  + 


H  r  X 


c  t  e 


H  P 


(C-4) 


This  can  be  reduced  as  follows 


*  c  o/T  0  T  c  ttW 

0  "  1  +  -EI -  -Q  (1+t])(o+1-1t)  +  -2-2 — 

HrXSp  H  r  XS 


(l+Tl)(a+l-iT) 


c  T  0  Tea*? 

<m) 


2  H  r  X 


(c-5) 


Substituting  the  definition  of  W  into  (C-5),  and  recognizing  that  —  =|  , 

rX  2 

end  simplifying  yields 


1  _  x  +  CpT3y  fo  [r)(a+l)  +  T  +  i  (ct+I-tjt)  ]  +  ^22 

a  '  /-STT  rt  n  M  H 


c  T  0„ 


2H  S  0 


H  3 


T  c  7  (T  -  T  )  r 


(c-6) 


=  1+02  [r)(o+l)  +t]  +  y  (a+i  -  |)  +  G 


+i  [G^  (o+1-t]t)  -  Yt] 


(C-7) 
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c  7T  0- 

where  3  =  PS  — 

wnere  ^  2  H  S  0 

2c  (T  -  T  ) 
Y  -  P  s  0 
1  “  H  S 


and 


G  = 


c  T  0- 
p  s  0 


H  0 


APPENDIX  D 


APPROXIMATION  FOR  0Q  -  THE  THERMODYNAMIC  LIMITED  CASE 


The  Clausius-Clapeyron  equation  is 


dp. 

dTp  RT  2 
s 

The  perturbed  pressure  and  surface  temperature  are 

T  =  T  [1  +  6  e(i+Tl)ait] 
s  s  o 

P-?  [1  +  *  e^1+I)  ^  J 

Substituting  (2.4)  into  (2.23)  yields 

( i+T)  )o>t 


P0«Ui  (i+Tj)  e 
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TABLE  5  (Continued) 


(a)  In  this  column,  the  abbreviation  "atm"  has  been  used  to  indicate 
that  the  burner  exhausted  to  the  atmosphere. 

(b)  In  this  column,  "C"  represents  chuffing;  "O"  represents  oscillatory 
burning;  "S"  represents  stable  burning;  ,'E"  represents  that  the 
sample  extinguished;  and  "U"  represents  uneven  or  irregular  burning. 

(c)  The  pressure  was  so  irregular  that  an  approximation  of  a  mean 
pressure  range  would  be  meaningless. 

(d)  L*-f requency  data  were  obtained  from  this  run. 

(e)  The  chamber  pressure  was  sufficiently  low  that  it  is  doubtful  that 
critical  flow  was  maintained  in  the  nozzle. 

(f)  A  record  of  this  run  was  not  obtained  with  the  tape  recorder. 
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irregular  burning. 

(c)  The  pressure  was  so  irregular  that  an  approximation  of  a  mean  pressure  range  would  be  meaningless. 

(d)  L*-frequency  data  were  obtained  from  this  run. 

(e)  A  record  of  this  run  was  not  obtained  with  the  tape  recorder. 


G  PROPELLANT  DATA  FROM  LOW-L* 


129  - 


JO 

CO 


co  d 

25 


42  C 

U 


a  d 


N 

N 

£ 


«  a 


44 

d 

(0 

r4  3 
8.8 
£ 


c| 
5  3 

25 


o 

CO 


0)  o 

0) 

o  h4 

o 

C  3 

d 

o  C 

o 

*3  T3 

-o 

4)  s  <u 

c-N  3  /-N 

44  TD  44 

•3  44  T3 

44  W  44 

W  44  W 

3  n  3 

3  i*s 

e  ^^6 

u  o  J:  u 

w  w  ©  w 

O  Pd 

l  l 

ouoo 


Pd  Pd 

l  I 

O  O  ©  O  U 


©  m  m 

r4  (M  H  H  M 


os  fl  B  as 

•  M  44  • 

41  3  <j* 


o  o  m  m  o 
miorsM^ 
•  •  •  •  • 
o  o  ©  ©  © 


CO  i-4  ©  ©  CM 
co  co  m  m  cm 


©  vo  O 

N  H  N  H 


SB  as  B  os 

44  •  44  • 

««*«** 


$££$$ 


m  m  in  ©  oo 

N  N  N  irin 


0) 

CO 

<98 

08 

08 

08 

£  8 

8 

i 

g 

>  -1 

5 

i4 

•3 

•H 

j>  ** 

44 

44 

44 

44 

44 

d 

3  »-4 

i-4 

CO 

CO 

H 

H 

»-4 

i-4 

0) 

(g 

m 

ttJ 

3 

co 

10 

If 

m 

•v  4 

u 

a 

B 

u 

u 

u 

u 

0 

<0  Q> 

Q J 

5 

i4 

3 

V 

3 

3 

i4 

O  > 

> 

44 

44 

U 

> 

> 

> 

44 

d  v 

3 

d 

3 

4) 

0) 

3 

o  <A 

«o 

<o 

n 

o 

08 

08 

08 

3 

*d  *o 

•3 

*3 

*3 

•3 

T3 

TJ 

•3 

*3 

*3 

0)  01 

/—s  0> 

0) 

01 

3 

0) 

4) 

3 

3 

3 

44  44 

T3  44 

44 

44 

CM 

44 

44 

44 

44 

44 

44  44 

W  44 

44 

44 

44 

44 

44 

44 

44 

44 

3  3 

^s  3 

3 

3 

3 

3 

3 

3 

3 

3 

66^66  S666 


Pd  p  p  © 
T  T  II 

o  a  o  ©  o 


VO  I  H  H  OS 
r4  CM  CM  CM 


Os  Os  B  O'  ^ 

•  •  44  •  • 

'J'  O  nt  N 


g  ^  vQ  Hi  H 


¥  ¥ 

©  ©  © 

I  I  I 

a  a  ©  © 


9999 

©  ©  ©  © 


44  44  44 


r-  «  O 
CM  ^  CM 


I  I 

I  CM  I 
CM 


oo  co  co 


OS  <f  OS  OS 
•  •  •  • 
N  ^  ^ 


n  n  n  <n 

iO  vfi  vO  vO 


r4  H  ©  H  H  H 


1-4  m  ©  m  m 
co  cm  oo  cm  CM 


as 


Os  m  Os  Os  Os 
•  •  •  •  • 
•a-  cm  <r  <t 


oo  co  co  ^ 

h.  rs  r*.  ro  co 


«-4  m  cm  m  cm 
CM  CM  CM  CM  CM 


g>  ©  ©  © 


o  _  _  _ 

m  <?S  m  m  m 

CM  CM  CM  CM  CM 


©  ©  ©  ©  © 

m  m  m  <n  m 

CM  CM  CM  CM  CM 


©  ©  ©  ©  © 

n  m  m  m  m 

CM  CM  CM  CM  CM 


©  ©  ©  © 
m  m  m  m 
cm  cm  cm  cm 


©  ©  ©  ©  © 

m  m  m  m  m 

CM  CM  CM  CM  CM 


<t  Mf 

OS  H  O'  OS  H 

I  I  I  I  I 

o  o  ©  ©  o 


0)  o 

m  vio  n  ^ 
H  ^4  »-4 

rJ  n  »4  N 
l  CM  I  I  CM 
O  I  ©  ©  » 

r-4  OS  H  i-4  Os 
co  sf  n  n  -s^ 


•O  N 

Os  O'  H  OS  H 
I  I  I  I  I 
©  ©  ©  ©  © 


f-4  *-4  OS  «-4  «-4 
till  » 
©  ©  ©  ©  © 


CM 

VO 

H 

i-4 

n 

3 

s-4  r* 
CM 

CO 

OS 

A 

m 

© 

i4 

vO 

i-4 

X-N 

3 

w 

§ 

CM 

3 

w 

Os 

i-4 

CM 

1 

-0-  r-l 

i-4 

»-4 

»-4 

H 

m 

m 

O 

1 

CM  1 

1 

1 

I 

1 

CM 

co 

CM 

CO 

o 

O 

1  © 

o 

o 

© 

o 

1 

i 

1 

OS 

t 

Os 

i-4 

i-4 

OS  r4 

^4 

i-4 

i-4 

1-4 

Os 

Os 

CO 

CO 

CO 

-3- 

<1- 

CO 

-3- 

<r 

■«0  ^  "O  sf 

f4  ^4 

I  I  I  I 

©  ©  ©  © 


? 


^  «jf 

H  H  H  i-4 
I  I  I  I 

ouoo 


i©  —>  <t 

in  m  m  <© 

??  a? 

o  ©  i  o 

^4  1-4  OS  ^4 

-tf  'S'  ^ 


49-254  G-14 


TABU!  8  (Continued) 


-  130  - 


•o 

V 

m 

3 

(A 

u 

S 

W 

o 

qj 

44 

44 

a 

u. 

•»4 

0) 

to 

r-4 

44 

o 

t3 

c 

CQ 

B 

U 

d 

o 

e 

t-4 

0) 

> 

XJ 

*o 

T3 

0) 

0) 

0> 

<y 

i-4 

■ 

<A 

LW 

44 

o 

C 

44 

44 

o 

A 

T3 

3 

0 

o 

0) 

/•s 

60 

V 

44 

6 

AN 

6 

w 

44 

0) 

✓"S 

3 

AS 

AN 

an 

/■N 

44 

6 

o 

w 

00 

w 

6 

6 

i4 

s*<* 

•i4 

i4 

NA 

•r4 

*H 

w 

•o 

CO 

3 

a 

» 

s 

u 

to 

* 

r-4 

O 

a  O 

o 

to 

Ed 

O 

CO 

© 

© 

C/3 

CO 

3 

w 

o 

o 

o 

D 

a 

u 

a 

© 

00 

CO 

CO 

o 

u 

o 

o 

u 

O 

O 

o 

C/3 

V 

Pd 

0) 

u 

0  ^ 

n  4)  m 

/-S  />  A  vX> 

B  M  rl 

44 

44 

44 

LM 

44 

CM 

0)  P  « 

£  a  £ 

o 

w 

1 

W 

I 

1 

W 

W 

I 

A* 

LA 

ON 

ON 

CO 

CM 

vO 

o 

vO 

A- 

NO 

O 

o 

NO 

CM 

1 

1 

1 

1 

i-4 

»-4 

r4 

r4 

CM 

CM 

CM 

CO 

CM 

CM 

CM  CO  <r 

i-4 

rs 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

<7> 

ON 

ON 

ON 

• 

a 

a 

a 

a 

a 

a  oo 

00 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

44 

44 

44 

■U 

44 

44 

• 

o  S* 

Nf  <t 

•d  *d 

•d 

*d 

*d 

•d- 

to 

eg 

to 

<0 

CO 

to 

ON 

ON 

Pd  w 

o 

■A 

lA 

LA 

o 

o 

o 

LA 

LA 

NO 

ON 

CO 

o 

LA 

lA 

o 

o 

LA 

LA 

LA 

o  C 

o 

CM 

CM 

CM 

LA 

LA 

LA 

r-» 

o 

»-4 

vO 

LA 

p*- 

r** 

LA 

LA 

P^ 

r- 

r** 

• 

• 

• 

• 

£ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

O 

o 

i-4 

r4 

r4 

O 

o 

o 

O 

O 

o 

o 

o 

la 

LA 

lA 

CM 

LA 

LA 

LA 

00 

LA 

LA 

LA 

LA 

*d 

MD 

H 

LA 

r-4 

r-4 

a  c 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

co 

CM 

CM 

• 

CM 

• 

CM 

d 

NO 

NO 

CO 

CO 

• 

CM 

• 

CO 

CO 

V 

o 

o 

O  O 

o 

O 

O 

00 

00 

00 

00 

00 

00 

oo 

°q 

00 

00 

OO 

00 

00 

LA 

LA 

lA 

LA 

LA 

LA 

LA 

CO 

CO 

co 

co 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

N  *»4  >i4 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

£  °w 

• 

• 

• 

• 

a 

•d' 

*d 

*d 

-d 

*d 

•d* 

r* 

r- 

P*- 

P*- 

P-. 

r— 

ON 

00 

<r 

•d 

-o 

H 

H 

H 

i-4 

*-4 

r4 

r-4 

r-4 

r-4 

ON 

ON 

OO 

ON 

i-4 

H 

8.8 

a 

o 

o 

o 

o 

u 

o 

o 

u 

o 

o 

o 

u> 

u 

o 

o 

u 

O 

o 

o 

O 

£ 

AN 

AN 

an 

A 

<u 

0> 

0) 

0) 

0) 

'w' 

w 

N-"' 

0) 

«o 

«d 

00 

i-4 

CM 

CM 

CM 

LA 

CM 

LA 

P- 

*d 

CO 

CO 

LA 

NO 

*d 

00 

r>» 

A- 

-4 

ON 

nO 

A* 

lA 

NO 

CO 

LA 

LA 

i-4 

8 

r^* 

c£ 

o 

o 

CM 

CM 

O 

CM 

i-4 

—4 

r4 

CM 

f4 

r4 

co 

r-4 

r4 

CO 

O 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

CM 

1 

1 

CM 

CO 

co 

D  0 

o 

o 

o 

o 

o 

O 

i-4 

O 

o 

O 

o 

O 

r-4 

o 

1 

^4 

H 

1 

1 

ON 

pH 

i-4 

r-4 

i-4 

^4 

r-4 

-4 

i-4 

r4 

r>4 

i-4 

i-4 

ON 

r4 

t-4 

ON 

ON 

£ 

-d- 

<r 

«a- 

-d 

-d 

*d 

•d- 

CO 

CO 

CO 

CO 

CO 

CO 

*d 

•d 

410-026  G-14  .228  .28  0.75  9.8  35-40  S 

410-027  G-14  .228  .28  0.75  9.8  40-42  S 

410-053  G-17  .228  .28  0.75  9.8  25-30  0-S 

410-07 1(e)  G-17  .228  .28  0.75  9.8  30  O-S 

410-054  G-17  .228  .23  0.75  9.8  30  S 


131 


T5  T> 

4)  a> 

CO  co 

00 

00 

GO 

GO 

GO 

3  3 

c 

c 

c 

C 

3 

•H 

•H 

u 

M 

u 

Xl 

l-l 

3 

TJ 

3 

-o 

3 

T3 

3 

T3 

3 

no 

0)  4) 

M  M 

a  cx 

co  co 

00 

GO 

GO 

GO 

GO 

4H 

3  3 

-0  *T3 

■u 

3  CO 

C  u 

C  CO 

C  co 

3  co 

MH 

C 

•H  c 

CO  o 

•*H  3 

CO  O 

‘H  C 

CO  O 

•H  C 
CO  O 

**H  3 

CO  o 

3 

X 

4)  4) 

•3  •a 

4)  a> 

pH  pH 

9 

«iH  *p4 

*iH  *H 
X*  4J 

•pH  **H 
Vh  -u 

*iH  hH 
U  4-» 

•pH  *H 
U  4J 

o 

CX  CX 

a 

CO 

0)  pH 

CO 

CL)  rH 

CO 

4)  pH 

CO 

Q)  pH 

CO 

0)  pH 

pH 

GO  GO 

O  0 

X  r-) 

Xl  pH 

Xl  pH 

U  rH 

X  pH 

s 

3  C 

3  *H 

CO  O 

3  *H 

CO  o 

3  -H 

CO  o 

3 

(0  O 

3  *H 

CO  O 

0 

CO 

cx  a. 
(X  £ 

U  >H 

4)  4) 

CO  CO 

CO  CO 

CO  CO 

CO  CO 

CO  CO 

^  *~\  0) 

a  3 

Mp*  w  O 

4)  O 
£ 

<U  0 

£ 

4)  o 
£ 

4J  0 
£ 

4)  0 
£ 

o  o 
u  n 

Q  Q 

C0  CQ 

9  9 

w 

0) 

D 

3 

CO 

pH 

CO 

c/a 

c/a 

C/a 

C/5 

O  CO  CO  CO  CO 

CO  CO 

3 

W 

o 

o 

o 

o 

to  o 

o  o  o  o  o 

O  O  CO 

gi  co  in  (A  n 

0) 

/—s 

P <5 

TJ 

<U 

C 

4J 

4) 

Xl 

5 

o 

co 

co 

o 

m 

o 

m 

co 

o 

O  to  o  o 

»n  pi  r-  •-< 

©  o 
mo© 

in  pH  i-H 

o  o  o 

ff  T 

c 

3 

w 

w  3  » 

U  (0  p- 

CM 

in 

CM 

o 

CO 

o 

CO 

vO  oo 

CO  CM 

O  O  O  CO  o 
cm  co  m  m  -a 

©  m  in 
m  <r  <* 

oommoo 
n  <t  <t  n  in 

00 

3 

IS 

/-^s  (0 

<Q  *H 

w0  a 

atm. 

8*6 

• 

00 

• 

O' 

7.4 

9.8 

B  B  oo  co  oo 

4J  4J  •  •  • 

co  co  as  as  as 

co  oo  oo 
•  •  ■ 
as  as  as 

oo  oo  oo  oo  oo 
•  •  •  •  ■ 
Os  Os  Os  Os  as 

*•> 
JC  c 

© 

m 

o 

m 

O 

m 

o 

m 

o  o 
in  m 

o  o  o  o  o 
m  m  m  in  m 

©  ©  © 
m  in  m 

o  ©  o  o  m 
m  m  m  m 
. 

O  *H 
•J  ^ 

O 

d 

o 

© 

©  o 

o  d  o  d  o 

o  ©  o 

o  o  ©  ©  © 

X“N 

cx  3 

H 

CO 

m 

CM 

CM 

CM 

os 

pH 

O'  vO 

^H  pH 

Ht  00  00  00  GO 

^  n  n  n  n 

.38 

.34 

.34 

pH  ph  oo  in  m 

CO  CO  CM  CM  CM 

w 

• 

4) 

*-H  •  /-s 

N  CO  C 
M  *H  -H 

Os 

O' 

O' 

O' 

as  as 

as  as  as  as  as 
o  o  o  o  o 

as  os  as 
©  ©  © 

as  as  as  as  as 

o  o  ©  o  © 

CM 

CM 

CM 

CM 

CM  CM 

CM  CM  CM  CM  CM 

CM  CM  CM 

CM  CM  CM  CM  CM 

o  a  ^ 
SB 

• 

* 

• 

a 

4-1 

c 

CO 

^  a> 

pH  *o 

Is* 

pH 

r-* 

H 

r-* 

H 

r- 

p4 

H  pH 

sr 

OS  OS  pH  H  H 

^  sf 

pH  H  H 

•tf  ^ 

pH  pH  pH  ^H  ^H 

8.3 

1 

o 

1 

o 

o 

o 

o  o 

o  o  o  o  o 

o  o  o 

o  o  o  o  u 

£ 

/-s 

/-N 

/■"S 

p— s 

4) 

0) 

0) 

4)  4) 

N-/ 

w 

w  w  4) 

m  pH  ^ 
m  CM  CM 
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«H  pH  O'  *"H 

o  ©  © 

»H  pH  ^H 

iiii© 
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'd- 

M  Cl  <J 

^  ^ 

<t  ~s-  ^  ^  ^ 

^ _ A 


410-0O5  G-17  .201  .44  0.50  atm.  —(f)  C-E  Chuffed  twice 

311-139  G-9  .201  .44  0.50  atm.  20-40  S 

410-072  G-17  .2")1  .41  0.50  atm.  60  S 
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TABLE  8  (Continued) 

(a)  In  this  column,  the  abbreviation  "atm"  has  been  used  to  indicate  that 
the  burner  exhausted  to  the  atmosphere. 

(b)  In  this  column,  "C"  represents  chuffing;  "0"  represents  oscillatory 
burning;  "S"  represents  stable  burning;  "E"  represents  that  the 
sample  extinguished;  and  "U"  represents  uneven  or  irregular 
burning. 

(c)  The  chamber  pressure  was  so  low  that  critical  flow  was  not  maintained 
in  the  nozzle. 

(d)  A  record  of  this  run  was  not  obtained  with  the  tape  recorder. 

(e)  L*-frequency  data  were  obtained  from  this  run. 

(f)  The  pressure  was  so  irregular  that  an  approximation  of  a  mean 
pressure  range  would  be  meaningless. 

(g)  A  thin  slab  of  F  propellant  approximately  1/16  inch  thick  was 
bonded  to  the  sample  in  an  effort  to  obtain  uniform  ignition. 

(h)  A  thin  slab  of  F  propellant  with  a  hole  in  the  center  was  bonded 
to  the  wall  of  the  burner  approximately  1/2  inch  from  the  sample 
for  ignition  purposes. 

(i)  The  minimum  pressure  of  the  pressure  oscillations  was  probably 
below  the  critical  pressure  needed  to  maintain  critical  flow 
in  the  nozzle. 


GB  PROPELLANT  DATA  FROM  LOW-L*  BURNER 
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TABLE  10 

CRITICAL  L*  AS  A  FUNCTION  OF  PRESSURE  -  F  PROPELLANT 


Run 

Number 

Pressure 

(psia) 

Cr itica 1 

L*  (in) 

42-063 

45 

9.3 

4/ -291 

34 

10.1 

410-161 

38 

15.0 

311-131 

48 

7.3 

47-251 

53 

7.4 

47-241 

53 

7.9 

47-252 

45 

9.5 

310-231 

39 

9.0 

410-146 

57 

10.0 

410-071 

50 

10.0(a) 

42-052 

80(b) 

4.4 

410-166 

88 

5.0 

410-163 

80 

5.8 

410-201 

56 

5.4 

311-133 

95(b) 

4.5(a) 

311-135 

85 

5.4 

(d)  Approximated  without  the  aid  of  the  planimeter. 

(b)  The  mean  pressure  at  which  the  propellant  burned  stably  after 
chuffing  once. 
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TABLE  11 

CRITICAL  L*  AS  A  FUNCTION  OF  PRESSURE  -  TF  PROPELLANT 


Run 

Number 

Pressure 

(psia) 

Critica 1 

L*  (in) 

48-261 

35 

20.9 

48-019 

48 

21 

47-282 

69 

15.4 

47 -254 

53 

18.8 

47-284 

50 

19.9 

51-042 

48 

20.8 

48-256 

68 

23(a) 

51-044 

130 

7.6 

48-011 

82 

12.0 

48-015 

110 

7.7 

48-016 

95 

6.4 

51-045 

115 

7.8 

48-017 

108 

7.8 

47-221 

122 

7.7 

48-264 

170 

6.5 

48-263 

140 

7.2 

(a)  Approximated  without  the  aid  of  the  planimeter 
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TABLE  12 

CRITICAL  L*  AS  A  FUNCTION  OF  PRESSURE  -  XF  PROPELLANT 


Run 

Number 

Pressure 

(psia) 

Critical 

L*  (in) 

48-266 

67 

21(a) 

49-027 

78 

9.9 

49-026 

95 

13(a) 

49-025 

75 

14.5 

48-269 

90 

10.5 

49-021 

120 

9.2 

49-022 

80 

12.6 

49-023 

170 

6.6 

49-024 

122 

8.4 

(a)  Approximated  without  the  aid  of  the  planimeter 
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TABLE  13 

CRITICAL  L*  AS  A  FUNCTION  OF  PRESSURE  -  G  PROPELLANT 


Run 

Number 

Pressure 

(psia) 

Critical 

L*(in)(a) 

39-233 

30 

24(b) 

310-163 

30 

22.7 

49-306 

35(c) 

25 

410-053 

32 

25 

410-071 

32 

23 

410-214 

30(c) 

12 

410-176 

34(d) 

14 

411-061 

42(d) 

15.3 

411-062 

34(d) 

16 

410-191 

34(d) 

17 

311-075 

30(e) 

7.1 

311-074 

36 

13.1 

410-022 

53 

8.4 

410-065 

60(f) 

3.5 

(a)  Determined  at  the  point  where  oscillations  cease. 

(b)  Approximated  without  the  planimeter. 

(c)  Pressure  uncertain  due  to  drift  in  transducer  output. 

(d)  The  mean  pressure  at  which  the  propellant  burned  stably  after 
oscillating. 
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TABLE  14 

L* -FREQUENCY  DATA  -  F  PROPELLANT 


Run  No. 

L*(in) 

Frequency 

(cps) 

Pressure 

(psia) 

1/S 

7 

49-181 

17  (a) 

50  (b) 

36 

42 

9.1 

19  (a) 

43  (b) 

35 

40 

8.2 

42-063 

7.2(a) 

121 

38 

95 

20.7 

9.0 

103 

44 

64 

14.9 

47-291 

10.1 

94  (b) 

34 

76 

18.2 

410-161 

15.0 

75 

38 

45 

12.8 

311-131 

7.3 

125 

,  48 

71 

16.4 

410-146 

8.1 

131 

57 

52 

13.9 

9.1 

115 

56 

48 

12.6 

10.0 

120 

53 

46 

13.8 

411-071 

10.0(a) 

105  (b) 

50 

50 

13.2 

42-052 

4.7(a) 

255 

85 

56 

17.1 

410-166 

5.0 

236 

88 

51 

15.3 

410-163 

4.2 

260 

70 

79 

21.8 

5.8 

226 

80 

49 

16.3 

311-135 

5.4 

241 

85 

49 

16.2 

410-147 

5.7 

250 

98 

39 

14.2 

410-211 

3.5(c) 

274 

80 

81 

19.8 

410-212 

5.0(c) 

225 

75 

61 

17.5 

(■3)  Approximate. 

(b)  Measured  from  less  than  four  complete  pressure  oscillations. 

(c)  Calculated  from  the  dimensions  of  the  remaining  propellant  after 
extinguishment . 
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TABLE  IS 

L* -FREQUENCY  DATA  FOR  TF  PROPELLANT 


Run  No. 

L*<in) 

Frequency 

(cps) 

Pressure 

(peia) 

1/5 

7 

n  (®) 

47 -286 

10  (b) 

91 

56 

52 

12.0 

•  •  • 

16  (b) 

63 

28  (c) 

62 

15.8 

•  •  • 

47 -287 

15  (b) 

74 

50 

38 

10.7 

0.041 

48-261 

15.7 

71 

37 

49 

13.8 

.062 

16.6 

67 

39 

44 

12.4 

.046 

17.4 

69 

39 

42 

12.8 

.050 

17.8 

71  <d) 

39 

41 

13.2 

•  •  • 

18.6 

62 

39 

40 

11.5 

.038 

20.9 

58 

35 

39 

11.9 

.000 

48-019 

17.6 

66 

51 

32 

9.4 

.040 

21 

65 

48 

29 

9.9 

.028 

47-282 

11.4 

91 

62 

42 

10.9 

.052 

12.3 

100 

73 

33 

10.3 

.073 

13.4 

93 

70 

31 

9.9 

.072 

15.4 

91 

69 

27 

9.8 

.035 

47 -284 

18.7 

66 

47 

33 

10.3 

.078 

19.3 

67 

52 

29 

9.4 

.088 

19.9 

67 

50 

29 

9.7 

.036 

51-042 

20.8 

74 

48 

29 

11.3 

•  ■  • 

48-256 

23  (b) 

72 

68 

19 

7.9 

.000 

51-044 

4.2 

194 

105 

68 

14.0 

•  •  • 

4.8 

214 

125 

51 

13.2 

•  •  • 

7.6 

209 

130 

31 

12.4 

•  •  • 

48-011 

8.4 

134 

70 

50 

14.3 

.135 

9.6 

133 

80 

38 

12.3 

.080 

12.0 

124 

80 

31 

11.5 

-.028 

14.8 

101 

58 

29 

11.0 

-.063 

51-043 

10.6 

128 

92 

31 

10.5 

•  •  • 

48-015 

4.1 

256 

130 

58 

15.2 

.106 

5.0 

214 

120 

51 

13.7 

.125 

6.0 

209 

120 

42 

13.4 

.082 

7.7 

190 

110 

36 

13.1 

.000 

51-045 

5.9 

161 

85 

59 

14.2 

•  •  • 

6.6 

178 

100 

46 

13.5 

•  •  • 

7.8 

189 

115 

34 

12.5 

•  •  • 

8.7 

196 

120 

29 

12.5 

•  •  • 

11.2 

168 

110 

22 

11.6 

•  •  • 
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TABLE  15  (Continued) 


Run  No. 

L*(in) 

Frequency 

(cps) 

Pressure 

(psia) 

1/5 

7 

n  (a) 

48-017 

7.8 

169 

108 

36 

11.9 

.111 

48-264 

3.4 

342 

215 

43 

12.7 

.145 

4.0 

316 

200 

39 

12.4 

.140 

4.7 

322 

200 

33 

12.6 

.135 

6.5 

237 

170 

28 

10.9 

.084 

48-263 

6.3 

224 

140 

35 

12.3 

.156 

7.2 

225 

140 

30 

12.4 

.138 

(a)  It 

was  not  poss 

ible  to  obtain 

approximations  for 

the  growth 

term 

on  all  of  the  runs. 


(b)  Approximate. 

(c)  The  minimum  pressure  of  the  pressure  oscillations  was  probably 
below  the  critical  pressure  needed  to  maintain  sonic  conditions 
in  the  nozzle. 

(d)  Measured  from  less  than  four  complete  pressure  oscillations. 


TABLE  16 


L* -FREQUENCY  DATA  FOR  XF  PROPELLANT 


Run  No. 

L*(in) 

Frequency 

(cps) 

Pressure 

(psia) 

1/6 

7 

48-266 

21  (a) 

65 

67 

20 

6.6 

49-027 

9.9 

109 

78 

38 

9.6 

49-026 

13  (a) 

89 

95 

25 

6.9 

49-025 

14.5 

77 

75 

37 

7.0 

48-269 

6.3 

178  (b) 

115 

45 

12.1 

6.8 

173 

112 

43 

12.0 

7.3 

173 

105 

42 

12.6 

7.9 

151 

100 

40 

11.5 

8.8 

143 

95 

37 

11.1 

9.8 

130 

92 

34 

10.3 

10.5 

134 

90 

32 

10.7 

49-021 

7.7 

133 

130 

34 

8.4 

9.2 

141 

120 

30 

9.3 

49-022 

9.7 

107 

85 

36 

8.8 

10.8 

109 

85 

32 

9.0 

12.6 

105 

80 

29 

9.1 

49-023 

3.6 

260  (b) 

195 

58 

12.6 

4.2 

236  (b) 

190 

50 

11.6 

5.2 

237 

180 

41 

12.1 

6.6 

200 

170 

33 

10.5 

49-024 

6.2 

197 

140 

40 

11.7 

7.7 

190 

125 

35 

12.3 

8.4 

181 

122 

33 

11.8 

(a)  Calculated  from  the  initial  L*. 

(b)  Measured  from  less  than  four  complete  pressure  oscillations 
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TABLE  17 

L* -FREQUENCY  DATA  FOR  G  PROP ELL/ 


Run  No.  L*  (in) 


Frequency  Pressure 

(cps)  (psia) 


49-231 

0.5(a) 

49-241 

10.5(a) 

49-242 

18  (a) 

410-194 

29  (a) 

49-308 

42  (a) 

49-308 

52  (a) 

410-052 

54  (a) 

410-064 

64  (a) 

410-073 

70  (a) 

311-153 

4.0(d) 

5.0 

6.5 

310-163 

4.5 

18.9 

311-154 

7  (d) 

8  (a) 

311-155 

9  (d) 

21  (d) 

49-306 

22 

410-071 

23.3 

410-214 

10.5(d) 

11.5 

410-176 

12  (d) 

13.2 

411-061 

14  (d) 

14.7 

411-062 

15  (d) 

16.0 

410-191 

16.7 

311-075  5.6 

7.1 
13.0 
8.0 


fflj 

|y| 

IqI 

fnj 

lyj 

|Q] 

311-074 

49-305 
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TABLE  17  (Continued) 

(a)  Calculated  from  the  dimensions  of  the  remaining  propellant  after 
extinguishment . 

(b)  Measured  from  less  than  four  complete  pressure  oscillations. 

(c)  Uneven  pressure. 

(d)  Approximated  from  the  pressure  trace. 

(e)  The  minimum  of  the  pressure  oscillations  was  probably  below  the 
pressure  needed  to  maintain  sonic  conditions  in  the  nozzle. 
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TABLE  18 

L* -FREQUENCY  DATA  -  GB  PROPELLANT 


Run  No. 

L* 

(in) 

Frequency 

(cps) 

Pressure 

(psia) 

1/5 

r 

53-066 

20 

(a) 

23 

27 

210 

24.4 

53-067 

30 

(b) 

22 

18 

223 

37.1 

31 

(b) 

30 

25 

148 

34.7 

36 

(b) 

25 

25 

131 

28.9 

53-068 

30 

(a) 

20 

23 

132 

25.4 

40 

(a) 

20  (c) 

23 

126 

25.4 

53-063 

12 

(a) 

55 

49 

179 

29.8 

53-064 

14 

(a) 

48 

46 

162 

27.5 

(a)  Calculated  from  the  dimensions  of  the  remaining  propellant  after 
extinguishment . 

(b)  Approximate. 

(c)  Calculated  from  less  than  four  complete  pressure  oscillations. 


APPENDIX  G 


THE  EFFECTS  OF  THERMAL 
RADIATION  ON  THE  ANALYSIS  OF  CHAPTER  II 

The  experimental  results  discussed  in  Chapter  IV  indicate  the 
possibility  of  a  different  mechanism  causing  NAI  in  G  propellant  than 
in  the  other  propellants.  It  vas  pointed  out  that  the  addition  of 
2  per  cent  carbon  black  to  the  composition  of  G  propellant  (to  make 
GB  propellant)  seemed  to  remedy  the  discrepancy  (i.e.,  data  from  GB 
propellant  firings  seemed  to  follow  the  same  trend  as  data  from  the 
other  propellants).  It  vas  concluded  that  at  least  part  of  the  reason 
for  the  discrepancy  was  due  to  penetration  of  thermal  radiation  into 
G  propellant.  Consequently,  an  attempt  was  made  to  include  Jhese 
effects  in  the  mathematical  analysis. 

Cantrell,  Me  _ure,  and  Hart  (7)  have  studied  the  effects  of 
thermal ^radiation  on  the  acoustic  response  of  solid  propellants. 
Although  they  were  interested  mainly  in  an  acoustic  phenomenon,  the 
qualitative  results  of  their  analysis  can  be  applied  to  the  problem 
at  hand,  llieir  analysis  indicated  that  radiation  Bhould  be  responsible 
for  a  considerable  increase  in  the  propellant  response  for  the  low 
frequency  regime  when  the  mean  burning  rate  of  the  propellant  is  rela¬ 
tively  low.  They  also  concluded  that  for  a  rather  typical  propellant, 
their  analysis  would  be  valid  for  burning  rates  down  to  approximately 
0.13  inches/second.  For  slower  burning  rates  the  effects  of  thermal 
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radiation  would  exceed  the  limits  of  their  assumptions .  The  burning 
rate  of  G  propellant  in  the  range  of  interest  is  less  than  0.10  inches/ 
second.  Therefore,  it  might  well  be  expected  that  rsuliation  should 
have  a  strong  effect  on  the  unstable  behavior  of  the  propellant. 

In  order  to  include  the  effects  of  thermal  radiation  it  was 
assumed  that  radiation  penetrating  the  solid  phase  would  obey  Beer's 
law  as  adapted  to  thermal  radiation  [11,  p*  383]*  This  can  be  written 


f 


r 


(G-l) 


where  f  is  the  radiant  flux  at  any  point  in  the  solid,  f  is  the 
r  rB 

radiant  flux  at  the  surface,  and  6  is  the  absorption  coefficient. 
Nomenclature  specific  to  this  Appendix  only,  has  been  included  at  the 
end  of  the  Appendix.  The  unsteady  heat  conduction  equation,  derived 
with  the  inclusion  of  this  term  to  account  for  radiation  effects,  is 


&T  ,  d2T  -5x 

P  c  xr  *  k  — n  +  P  rc  x*  +  6f  e 
P  P  «t  gx2  p  p  ox  rs 

Solving  the  steady-state  form  of  this  equation  yields 


(0-2) 


T  -  TQ  + 


(*. '  ♦  1$)  •**  -  & 


(0-3) 


where  f  » 


,  and  6^-7^- 
T  rpc  r 


Differentiating  (0-3)  gives  the  thermal  gradient  in  the  solid  including 
the  effect  of  radiation  as 


-  ii)8  - 


dT 

dx 


_  /„  flT8\  -2Xx  6flT8  -&x 

2X  Ts  -  To  -  5^1  6  +  VT  6 


(0-4) 


The  perturbation  of  the  incident  radiative  flux  at  the  surface  of 
the  propellant  is  defined  as  the  other  perturbed  quantities  have  been 
defined 

f  -  f  [1  +  4>  e(i+Tl)ait]  (a-5) 

rs  rs 

where  it  is  recognized  that  the  radiative  flux,  f  is  a  fraction  of 
the  total  incident  flux  at  the  surf  Me,  f.  Substituting  the  perturbed 
radiant  flux,  burning  rate,  and  temperature  (see  Equations  (G-5), 
(2.4a),  and  (2.4b))  into  Equation  (G-2),  subtracting  the  steady-state 
terms,  and  neglecting  second  order  perturbations,  yields  t..e  following 
second  order  differential  equation  in  terms  of  the  perturbed  quantities 


4f  +  =*  g  -  (LI)  |  ^  « 

dX  Cz  m  ax  rvrv/%  O’ 


Bf 


-5x 


ape  T 
P  8 


Substituting  Equation  (G-4)  into  (G-6)  results  in 


(0-6) 


+2xg  -  (ti+i)2  9  b  4X2p 
dx 


T  -  Ta  f. 

s  0  1 


-2Xx 


V1 


erf 


I  [gg—  + «]  e_6x 


Applying  the  boundary  conditions  that 

at  x  ■  0  9(0)  ■  Qr 


(0-7) 


and  at  x  ■  « 


0(oo)  m  0 
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to  Equation  (G-7)  and  solving  as  in  Chapter  II,  where  radiation  effects 
were  not  included,  gives  the  solution 


-  £>  -  v  c2>  .  Ci 


e'2Xx  +  C„  .-8x 
c. 


(0-8) 


where  C. 


s 

7(n+i) 


c 


2 


-  fi  (57^1*1) 

Bl'1  '  T&1  (T1  +  15 


Equation  (G-8)  is  analogous  to  Equation  (2.1l)  of  the  main  text. 

The  incentive  for  deriving  an  equation  such  as  (G-8)  was  to  permit 
the  integration  of  the  integral  in  Equation  (2.6).  Recalling  that 
the  radiant  flux  is  a  fraction  of  the  total  energy  flux,  f,  it  can 
then  be  seen  that  Equation  (2.6)  will  be  the  same  even  when  the  effects 
of  radiation  are  considered.  For  convenience  this  equation  has  been 
repeated 


1 

3 


x  +  fo  +  CP  ^  ^5. 


(2.6) 


Utilizing  the  results  of  Equation  (G-8)  the  indicated  integration  can 
be  carried  out  as  in  Appendix  C.  The  results,  after  algebraic  manipu' 
lation,  are 


|  -  1  +  G  +  0,  h(o+l)  +  t+  i(o+l-i)T)]  -  T,[S  -  2(o+l-it)]  (G-9) 

0  X 


-  Y 


f  0 
1  1 


k  (?2-n)2  +  1 


(oT1  -t)  +  0  +  1  -  g|-  (r)2+l) 


'2  '26 


+  i  [r2  (gf-  +  ir0  -  T(Ti  +1)  *  tC0-*-1)] 
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vhere  Y_ 


:  i  TT  -T-  f,  "I 
JE_»  .s,-..0  +  __1_ 

b  L  j  6r1-! 


2c  * 
HS 


7cT 
P  5 
2HS 


72  - 


_i_  +  £ 

a^i  P 


Ifc^-l) 


Separating  Equation  (G-9)  Into  lta  real  and  Imaginary  parts  yields 


|  -  »•  <g>  * 1  r“  <|> 

f •%  * 


(G-10) 


where  Re  (|)  fi  1  +  G  +  h(o+l)  +  t)  -  Y3  [S  -  2  (o+l)] 


-  ^  *2 1»,  (g,  -  -o  * .  *  1  -  ^ 


In  (|)  ^  Gx  (a+l-T)t)  -  2  Y3t 


®2  [72  (2B_^  +1t) 


-  t(t!  +l)-T](o+l)  ] 


and  1r 


t  ♦ 
1  1 


(r2-Tir  +  1 


This  equation  is  analogous  to  Equation  (2.12)  and  can  now  he  combined 
with  the  results  of  the  perturbed  mass  balance  (Equation  (2.17)) 
obtaining  an  equation  analogous  to  Equation  (2.19)- 

This  was  done  and  the  results  were  programed  and  run  on  the 
IBM  70UO  computer.  However,  positive  results,  agreeing  with  the 
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experimental  data,  were  not  obtained  in  the  initial  effort, 
the  project  was  not  pursued. 


Therefore, 


NOMENCLATURE  FOR  APPENDIX  G 


ROMAN  LETTERS 


-k 


(*."T0)  fl 

f  +  ^l"1 


TT^iT 


-  fi  [  »j-i 


+  8 


V1  •  i«r  (n+1) 


radiant  flux 


re 


Tp  r  c 

8  p 


2c  T 

_ E.JL 

HS 


T8-T0  ,  fl 

ft  Bl'1 
>-  8 


GREEK  LETTERS 


7CT 
.  J2-S 
2HS 


^(8^1) 


6 

8, 


adsorption  coefficient 
06 


radiant  flux  perturbation  at  the 
surface 


dimensionless 


dimensionless 


m/e3 


dimens ionlesB 


dlmenslonles  s 


dimensionless 


dimensionless 

1/1* 

dimensionless 


dimensionless 
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GREEK  LETTERS 
(Cont'd) 


$> 


1 


it 


2 


+  i 
a^i  e 


(72-ri  )2+  1 


dimensionless 


dimensionless 


PAGKS  I54  THROUGH  158  OMITTED 
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TABLE  19 


Strand 


Run 

No. 

Propellant 

Type 

Length 

(in.) 

0 

(°c) 

P 

(psia) 

Rate 

(in. /see. ) 

3 

F-40 

2.0 

25 

45 

0.216 

5 

F-40 

2.0 

19 

34 

0.193 

9 

F-40 

2.0 

20 

53 

0.241 

12 

F-40 

2.0 

21 

32 

0.169 

13 

F-40 

2.0 

24 

28 

0.164 

15 

F-40 

2.0 

25 

38 

0.216 

17 

F-40 

1.98 

25 

43 

0.208 

21 

F-40 

2.0 

29 

39 

0.200 

22 

F-47 

2.0 

28 

18 

0.128 

23 

F-47 

2.0 

28 

22 

0.140 

24 

F-47 

2.0 

28 

28 

0.166 

25 

F-47 

2.0 

28 

33 

0.176 

26 

F-47 

2.0 

26 

37 

0.190 

27 

F-47 

2.0 

26 

43 

0.206 

28 

F-47 

2.0 

26 

47 

0.225 

29 

F-47 

2.0 

26 

51 

0.234 

30 

F-47 

2.02 

28 

57 

0.254 

31 

F-47 

2.02 

28 

64 

0.268 

32 

F-47 

2.0 

29 

67 

0.286 

33 

F-47 

2.02 

26 

66 

0.284 

34 

F-47 

2.03 

29 

14 

0.120 

35 

F-47 

2.00 

30 

15 

0.124 

36 

F-47 

2.02 

30 

17 

0.126 

38 

F-47 

2.00 

30 

17 

0.128 

39 

F-47 

2.00 

30 

20 

0.135 

4l 

F-47 

1.98 

31 

26 

0.160 

42 

G-12 

2.00 

27 

26 

0.0798 

43 

0-12 

2.02 

27 

17 

0.0642 

45 

0-12 

2.01 

29 

18 

0.0640 

46 

0-12 

2.02 

30 

21 

0.0720 

47 

G-12 

2.01 

31 

22 

0.0720 

48 

G-12 

2.00 

32 

23 

0.0737 

49 

G-12 

2.00 

32 

25 

0.0793 

50 

G-12 

1-99 

30 

26 

0.0790 

51 

G-12 

2.00 

30 

40 

0.0943 
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TABLE  9  (continued) 


Run 

No. 

Propellant 

Type 

Strand 

Length 

(in.) 

T0 

(#C) 

P 

(psia) 

Rate 
(in. /sec 

52 

G-12 

2.00 

29 

28 

0.0839 

53 

G-12 

2.00 

30 

30 

0.0841 

5^ 

G-12 

2.00 

31 

32 

0.0913 

55 

G-12 

2.00 

31 

34 

0.0876 

56 

G-12 

2.00 

28 

14 

0.0545 

57 

G-12 

1-99 

28 

15 

0.0592 

58 

G-12 

2.00 

29 

16 

0.0593 

59 

G-12 

2.02 

29 

17 

0.0615 

60 

G-12 

2.00 

31 

19 

0.0674 

61 

G-12 

2.00 

31 

36 

0.0887 

62 

G-12 

2.01 

32 

37 

0.0900 

63 

G-12 

2.01 

32 

4l 

0.0949 

64 

G-12 

2.01 

33 

48 

0.1091 

65 

G-12 

2.00 

30 

51 

O.IO83 

66 

G-12 

2.00 

31 

66 

0.119 

67 

G-12 

2.00 

31 

56 

0.105 

68 

G-12 

2.00 

31 

62 

0.112 

69 

G-12 

2.00 

31 

66 

0.116 

71 

G-12 

2.00 

32 

73 

0.129 

72 

F-69 

2.00 

32 

16 

0.127 

73 

F-69 

2.00 

29 

74 

0.302 

74 

F-69 

2.01 

29 

83 

0.335 

76 

F-69 

2.00 

30 

103 

0.373 

77 

F-69 

2.01 

31 

115 

0.375 

78 

F-69 

1.99 

31 

125 

0.396 

79 

F-69 

2.01 

31 

142 

0.411 

80 

F-69 

2.00 

30 

167 

0.481 

81 

F-69 

2.01 

30 

188 

0.501 

82 

F-69 

2.01 

30 

201 

0.506 

83 

F-69 

2.00 

30 

2l4 

0.563 

84 

TF-1 

2.00 

32 

29 

0.175 

85 

TF-1 

2.01 

32 

48 

0.217 

86 

TF-1 

1-99 

32 

72 

0.286 

87 

TF-1 

2.01 

32 

97 

0.344 

88 

F-69 

1.98 

29 

75 

0.292 
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TABLE  9  (continued) 


Run 

No. 

Propellant 

Type 

Strand 
Length 
(in. ) 

89 

F-69 

1.98 

90 

F-69 

1.98 

91 

F-69 

1.97 

92 

F-69 

2.00 

93 

F-69 

2.00 

94 

F-69 

2.00 

95 

F-69 

2.01 

96 

TF-1 

2.00 

97 

TF-1 

2.00 

98 

TF-1 

2.00 

99 

TF-1 

2.00 

100 

TF-1 

2.00 

101 

TF-1 

2.00 

102 

TF-1 

2.00 

103 

TF-1 

2.01 

104 

F-69 

1-97 

105 

f-69 

1-97 

106 

f-69 

1.98 

107 

F-69 

1.98 

108 

TF-1 

2.00 

109 

TF-1 

1-99 

no 

TF-1 

1-99 

112 

G-12 

2.00 

113 

G-12 

2.01 

114 

G-12 

2.01 

115 

G-12 

2.00 

ll6 

G-12 

2.01 

117 

G-12 

2.01 

118 

G-12 

2.00 

119 

G-12 

2.00 

120 

G-12 

2.01 

121 

G-12 

2.00 

1’ 

F-69 

2.00 

2' 

f-69 

2.01 

3' 

F-69 

2.01 

iO  P  Rate 

(°C)  (psia)  (in. /sec.) 


30 

91 

0.338 

30 

105 

0.343 

30 

129 

0.386 

31 

162 

0.4S4 

31 

186 

0.476 

32 

229 

0.499 

32 

254 

0.523 

33 

15.3 

0.151 

34 

21 

0.148 

34 

39 

0.200 

34 

60.5 

0.247 

34 

88 

0.304 

35 

121 

0.348 

35 

145 

0.385 

35 

207 

0.439 

28 

14 

0.1145 

29 

15 

0.124 

29 

16 

0.126 

30 

17 

0.134 

34 

16 

0.139 

34 

21 

0.163 

34 

23 

0.176 

32 

94 

0.166 

32 

114 

0.186 

32 

136 

0.210 

32 

152 

0.208 

31 

69 

0.l4l 

31 

55 

0.128 

32 

40.5 

0.114 

33 

28.5 

0.0911 

30 

14 

0.0602 

30 

17 

0.0725 

33 

19 

0.167 

33 

15 

0.139 

34 

44 

0.245 
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TABLE  9  (continued) 


Strand 


Run 

No. 

Propellant 

Type 

Length 

(in.) 

0 

(°C) 

P 

(psia) 

Rate 

(in. /sec .  ) 

4' 

F-69 

2.01 

34 

80 

0.362 

6' 

F-69 

2.01 

35 

28 

0.195 

7' 

F-69 

2.00 

35 

94 

0.420 

8' 

F-69 

2.01 

35 

126 

0.455 

9' 

F-69 

2.00 

34 

69 

0.339 

10' 

F-69 

2.01 

34 

59 

0.302 

122 

G-12 

2.00 

30 

22 

0.0817 

123 

G-12 

2.01 

30 

20 

0.0753 

124 

G-13 

2.01 

34 

20.5 

0.0788 

125 

0-13 

2.01 

34 

14 

0.0612 

126 

G-13 

2.01 

34 

18 

0.0815 

127 

G-13 

2.00 

34 

43 

0.117 

128 

G-13 

1  -99 

32 

27 

0.100 

129 

G-13 

2.00 

33 

63 

0.164 

130 

G-13 

2.00 

34 

83 

0.181 

131 

G-13 

2.00 

33 

102 

0.192 

132 

G-13 

2.00 

33 

129 

0.216 

133 

G-13 

2.00 

33 

165 

0.242 

144 

G-13 

2.01 

29 

195 

0.251 

145 

G-13 

2.01 

30 

226 

0.260 

146 

G-13 

2.00 

30 

268 

0.266 

147 

G-13 

2.00 

32 

299 

0.273 

148 

G-13 

2.02 

31 

362 

0.296 

149 

G-13 

2.02 

32 

410 

0.427 

150 

G-13 

2.00 

28 

73 

0.151 

152 

G-13 

2.02 

29 

145 

0.205 

153 

G-13 

2.02 

29 

192 

0.267 

154 

G-13 

2.01 

30 

225 

0.258 

155 

G-13 

2.01 

30 

310 

0.295 

156 

G-13 

2.00 

31 

412 

0.325 

158 

F-69 

2.02 

27 

314 

0.587 

159 

F-69 

2.00 

27 

416 

O.638 

160 

F-69 

2.00 

2  8 

497 

0.669 

161 

TF-1 

2.01 

31 

309 

0-535 

162 

TF-1 

2.02 

32 

414 

0.601 
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TABLE  9  (continued) 


Run 

No. 

Propellant 

Type 

Strand 
Length 
(in.  ) 

T0 

(#C) 

P 

(psia) 

Rate 
(in. /sec 

l64 

XF-2 

2.01 

27 

13-5 

0.142 

165 

XF-2 

2.00 

28 

18 

2.153 

166 

XF-2 

2.01 

28 

29 

0.193 

168 

XF-2 

1.98 

26 

39-5 

0.224 

169 

XF-2 

1.98 

28 

54 

0.255 

170 

XF-2 

2.01 

24 

69-5 

0.279 

171 

XF-2 

1-97 

24 

90 

0.296 

172 

XF-2 

2.01 

24 

114.5 

0.358 

173 

XF-2 

1.97 

20 

143 

0.371 

174 

XF-2 

2.00 

20 

170 

0.385 

175 

XF-2 

2.00 

19 

15 

0.134 

176 

XF-2 

2.00 

20 

23 

0.146 

177 

XF-2 

2.01 

20 

221 

0.431 

179 

XF-2 

2.01 

23 

33 

0.199 

180 

XF-2 

2.01 

24 

49.5 

0.239 

182 

XF-2 

1.99 

24 

290 

0.453 

183 

XF-2 

2.01 

25 

383 

O.5O8 

184 

XF-2 

1.99 

25 

310 

O.561 

185 

XF-2 

1.99 

24 

13 

0.130 

186 

XF-2 

2.01 

25 

19 

O.161 

187 

XF-2 

2.00 

25 

27 

0.182 

190 

U-l 

2.01 

25 

13 

0.148 

194 

U-l 

2.51 

26 

19 

0.158 

195 

U-l 

2.50 

26 

23 

0.185 

197 

U-l 

2.49 

24 

38 

0.227 

198 

U-l 

2.52 

24 

53 

0.274 

199 

U-l 

2.51 

25 

74 

0.341 

200 

U-l 

2.51 

25 

104 

0.418 

201 

U-l 

2.51 

25 

156 

O.5O6 

202 

U-l 

2.50 

21 

13 

0.127 

204 

U-l 

2.52 

22 

13 

0.127 

206 

U-l 

2.53 

23 

15 

0.140 

207 

U-l 

2.51 

23 

232 

0.612 

I  f 
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13  ABSTRACT 

Non-acoustic  combustion  instability  has  been  examined  in  an  uncatalyzed,  a 
catalyzed,  and  two  aluminized  composite  propellants.  These  propellants  were 
studied,  burning  cigarette  fashion,  in  a  burner  capable  of  operating  at  values  of 
Lx  as  small  as  5  cm.  It  was  observed  that  the  rrequency  of  the  pressure  oscilla¬ 
tions  varied  with  the  value  of  L*,  frequency  decreasing  with  increasing  L*.  The 
data  were  correlated  by  plotting  frequency  versus  the  reciprocal  of  L*,  yielding 
a  series  of  constant  pressure  curves.  This  pressure  effect  was  eliminated  by 
using  dimensionless  variables,  allowing  all  of  the  data  for  a  given  propellant 
to  be  correlated  along  the  same  line. 

A  one-dimensional  model  is  proposed  that  considers  sinusoidal  perturbations, 
allowing  for  growth  ot  tne  aist-rcance.  The  pressure,  the  burning  rate,  the 
attributed  temperature  rr.  tne  propellant,  and  the  energy  flux  fro:.,  tne  burning 
gases  are  the  quantities  parturoed. 
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